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quaere  film  dampers  are  widely  used  components  in  modern  turbomachinery,  but 
their  behavior  is  not  well  understood.  This  program  was  Initiated  In  order  to 
extend  the  squeeze  film  damper  analytical  basis,  confirm  predictions  by  means 
of  a  test  program,  transfer  the  analytical  technology  to  a  gas  turbine  engine 
manufacturer,  integrate  it  with  existing  rotordynamic  tools  sad  develop  and 
demonstrate  a  damper  design  system. 
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SECTION  I 


INTRODUCTION 

With  trends  towards  higher  thrust- to-weight  ratios,  improved  fuel  consump¬ 
tion  and  increased  durability  in  advanced  gas  turbine  engines,  conflicting 
design  constraints  have  arisen.  Paragraph  3.3.8.10  of  MIL-E-5007D  states 
"the  engine  shall  be  free  of  destructive  vibration  at  all  engine  speeds 
and  thrust,  including  steady-state  and  transient  operation  throughout  the 
complete  operating  range  of  the  engine".  Meeting  these  needs  while  design¬ 
ing  lighter  and  more  flexible  rotors  and  cases  has  significantly  increased 
the  requirement  for  analytical  design  capabilities  that  provide  for  the 
precise  control  of  the  rotor  system  dynamics.  The  overall  intent  of  this 
program  is  to  provide  a  better  understanding  of  the  physical  phenomena  that 
affect  the  design  of  Squeeze  Film  Dampers  for  high  speed  rotating  machinery 
and  to  provide  an  Integrated  rotor/ damper  design  system  methodology.  The 
integration  of  the  squeeze  film  damper  and  rotor  dynamics  design  tools  is 
seen  as  being  crucial  to  the  development  of  rotor  systems  characterized  by 
predictable  and  reliable  operation. 

The  team  of  Mechanical  Technology,  Inc.  (MTI)  and  Detroit  Diesel  Allison 
(DDA)  is  addressing  this  task  by  means  of  a  two  phase  program: 

Phase  I  -  Large  Squeeze  Film  Damper  Technology 

Phase  II  -  Small  Squeeze  Film  Damper  Technology  and 
Cavitation  Effects 

This  first  interim  report  documents  the  Phase  I  analysis  and  design  systems 
developed  under  this  program  by  both  MTI  and  DDA.  At  the  present.  Phase  I 
analysis  tasks  have  been  completed.  Verification  testing  through  the  use 
of  a  high  speed  controlled  orbit  rig  is  in  its  final  stages.  Work  on  Phase 
II  has  just  begun,  in  parallel  with  the  tasks  in  the  latter  half  of  Phase  I. 

The  following  paragraphs  outline  the  contents  of  the  two  phases,  indicate 
where  the  program  currently  stands  and  what  will  be  reported  on  in  this  report. 

Phase  I,  Large  Squeeze  Film  Damper  Technology,  consists  of  nine  tasks.  During 
the  first  three  tasks  of  Phase  I,  MTI  is  to  develop  a  squeeze-film  damper 
component  analysis  tool,  and  then  verify  this  capability  through  high  speed 


testing  of  up  to  five  different  design  geometries  and  sealing  conditions. 
This  verified  component  technology  is  then  to  be  assembled  and  documented 
in  a  form  that  could  be  directly  transferred  to  the  design  systems  of  mili- 
tary  gas  turbine  engine  manufacturers .  To  ensure  this  transferability,  DDA 
is  to  accept  the  verified,  documented  technology  and  integrate  it  with  their 
rotordynamic  design  and  analysis  capability.  DDA  is  to  then  develop  inte¬ 
grated  design  procedures  which  will  be  applied  both  in  application  rig 
design  and  in  design  improvement  of  a  current  inventory  Air  Force  engine 
(the  TF41)  which  .powers  the  single-engine  LTV  A7D  Corsair  2.  KTI,  in 
parallel  with  this  technology  transfer  and  application,  will  oversee  the 
design  system  development  and  provide  program  management  and  technology 
documentation . 

Component  analysis  and  component  rig  design  proceeded  In  parallel.  Initial 
transfer  of  software  to  DDA  took  place  once  the  component  analysis  was  com¬ 
plete.  Component  tests  are  to  be  performed  and  comprehensive  design  charts 
developed. 

Comparisons  between  predictions  and  measurements  are  to  be  made  once  tests 
have  been  completed  on  the  controlled  orbit  rig  and  the  free  orbit  rig. 

Under  Task  IV,  any  needed  refinements  to  test  or  analysis  are  to  be  planned 
and  full  component  technology  transfer  to  DDA  can  be  completed  and  inte¬ 
grated  fully  with  the  DDA  rotordynand.es  analysis  capability.  In  parallel 
with  the  component  technology  transfer,  MTI  is  to  do  further  work  in  consoli¬ 
dating  the  component  analysis  software  into  the  component  design  system  for 
use  in  association  with  the  integrated  rotor/ damper  system  design. 

Utilizing  the  fully  developed  and  verified  damper  component  and  rotor  damper 
design  system  capability  and  DDA's  experience  in  military  gas  turbine 
design,  DDA  is  to  develop,  document  and  illustrate  the  procedures  which 
should  be  followed  to  most  effectively  apply  the  damper  design  technology. 
Work  will  then  start  on  a  design  manual.  In  parallel  with  this  reporting 
effort,  the  design  criteria  will  be  applied  to  a  current  Air  Force  inventory 
engine  (TF41) .  One  or  more  dampers  for  this  engine  will  be  designed  and 
the  benefits  of  these  design  modifications  will  be  analytically  demonstrated 


Application  rig  work  is  also  planned  as  a  verification  of  the  design  pro¬ 
cedures.  Tests  are  to  be  performed  on  a  blade-loss  test  rig,  which  will 
demonstrate  dampers  designed  both  for  normal  and  for  abusive  unbalance. 
Comparisons  of  test  results  and  predictions  will  be  made.  Evaluation  of 
predictive  capability  will  be  performed  as  well  as  evaluation  of  the  design 
procedures  used  in  developing  the  damper  components. 
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Phase  II,  Small  Squeeze  Film  Damper  Technology  and  Cavitation  Effects, 
consists  of  six  tasks  which  complement  the  Phase  I  effort.  The  investiga¬ 
tions  proposed  in  Phase  II  will  lead  to  further  state-of-the-art  squeeze- 
film  damper  advancement,  both  in  fundamentals  and  in  applications.  The 
main  thrust  of  the  work  planned  is  to  better  understand  the  phenomenon 
of  cavitation  in  squeeze-film  dampers,  to  conduct  parameter  studies  of 
various  damper  phenomena  using  existing  software  and  to  assess  the  influence 
of  small  engine  size  (in  the  500-  to  8500-lb  thrust  class)  on  engine  and 
damper  dynamics. 

The  first  three  tasks  of  Phase  II  are  devoted  to  a  detailed  study  of  cavi¬ 
tation.  The  present  analysis,  while  it  considers  cavitation,  does  not 
provide  for  an  in-depth  investigation  of  this  phenomenon  which  can  strongly 
influence  damper  behavior.  The  scope  described  in  this  effort  will  develop 
understanding  and  a  means  of  handling  cavitation  effectively  in  design 
analysis . 

The  fourth  task  involves  a  series  of  parameter  studies  which  have  been 
selected  to  allow  the  full  benefit  of  the  software  capabilities  developed 
under  the  present  Air  Force  program  to  be  exercised  and  exploited,  particu¬ 
larly  in  the  nonlinear  area. 

The  investigation  of  sizing  effects  in  the  fifth  task  is  important  because 
the  present  Air  Force  program  is  directed  at  larger  engines.  Although  the 
basic  technology  should  be  extrapolatable  to  smaller  engines,  careful  con¬ 
sideration  of  sizing  effects,  both  from  damper  considerations  and  from 
engine  dynamics  considerations,  is  to  be  carefully  assessed. 
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The  overall  program  Is  outlined  In  the  block  diagram  of  Figure  1-1, 
arranged  by  type  of  effort  rather  than  by  task.  The  heavily  outlined  blocks, 
damper  analysis  and  design  system  development  are  completed  whereas  the 
boxes  outlined  with  broken  lines, 

•  steady  state  verification 

•  transient  verification 

•  design  methods  application 

•  cavitation  analysis 

•  visual  observation  of  cavitation 

are  currently  in  progress,  some  being  nearly  completed. 

This  first  interim  report  covers  the  two  completed  portions  of  the  program 
effort.  The  squeeze  film  damper  film  analysis  is  given  in  detail  including 
the  boundary  conditions  assumed ,  solution  methods  the  computer  program 
system  and  example  calculations.  The  design  system  development  is  also 
covered  in  detail,  including  damper  analysis/rotordynamic s  tools  and 
integration  philosophy  and  the  interfacing  of  these  tools  into  a  system 
for  linear  and  nonlinear  calculations. 


Figure  1-1  Block  Diagram  of  Squeeze  Film  Damper  Program  Effort 


SECTION  II 


SQUEEZE  FILM  DAMPER  ANALYSIS 

This  Section  describes  analytical  basis  behind  the  squeeze  film  damper 
analysis  computer  program. 

In  its  most  general  form,  the  computer  program  calculates  response  of  the 
film  to  an  imposed  velocity  at  a  particular  location.  Circular  orbit  cal¬ 
culations,  where  the  velocity  direction  is  always  at  90°  to  the  eccentricity 
vector,  are  a  special  case  of  this  generalized  calculation.  A  finite 
difference  approximation  to  the  Reynolds  equation  is  used  which  provides 
algebraic  difference  equations  linking  pressure  values  at  each  point  of  a 
rectangular  mesh  on  the.  developed  damper  surface.  A  very  mature  state  of 
the  art  exists  for  this  portion  of  the  problem.  Film  rupture  is  handled  by 
setting  computer  subambient  pressures  to  ambient  (with  an  option  for  ambient 
pressure  to  be  set  at  a  nonzero  value) .  It  is  noted  that  the  choice  was 
made  to  treat  the  zero  pressure  gradient  boundary  condition  as  a  refinement, 
which  could  be  added  later  if  the  need  arose. 


The  specific  effects  which  are  included  in  the  computer  program  are: 

1.  The  finite  length  of  the  damper,  which  is  handled  by  means 
of  the  two-dimensional  finite  difference  solution. 

2.  Finite  orbit  size  of  eccentricity  is  an  implicit  part  of 
the  finite  length  Reynolds  equation  solution. 

3.  Finite  velocity  is  accounted  for  in  terms  of  its  effect 
on  cavitation.  In  combination  with  the  supply  pressure, 
the  velocity  determines  the  location  of  points  in  the  film 
where  the  pressure  changes  from  positive  to  negative.  By 
tracking  these  points,  finite  velocity  effects  are  accounted 
for.  It  is  noted  that,  in  the  special  case  of  zero  supply 
pressure,  the  location  of  the  boundary  conditions  is  defined 
strictly  by  the  velocity  direction;  therefore,  in  this 
special  (though  meaningful)  case,  fluid  film  forces  are 
directly  proportional  to  velocity  magnitude;  a  fact  to  be 
remembered  when  making  effective  use  of  this  computer  program. 

4.  Damper  clearance  and  diameter  are  implicit  in  the  physics  of 
the  problem  represented  by  the  Reynolds  equation. 

5.  Fluid  viscous  effects  are  also  Implicit  in  the  Reynolds 
equation. 
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Fluid  inertia  effects  are  accounted  for  by  a  perturbation 
solution  based  upon  the  Reynolds  number  (X) .  The  approach 
used  is  that  described  by  Lund  in  Reference  1.  The  result 
of  the  fluid  inertia  effects  calculation  is  a  linear  correc¬ 
tion  to  the  predicted  pressure  and  viscous  damping  results. 

Fluid  film  rupture  is  handled  normally  by  the  ir  or  half¬ 
film  type  approach  where  any  computed  subambient  pressures 
are  set  to  and) lent  pressure  and  ambient  pressure  can  be 
set  as  an  Input  parameter. 

Fluid  supply  conditions  are  treated  by  specifying  pressure 
to  be  equal  to  inlet  pressure  over  the  region  of  the  supply 
groove  or  hole.  Flow  boundary  conditions  and  inlet  and 
outlet  will  be  handled  by  satisfying  either  the  local  or 
global  relationships  between  flow  and  pressure  drop.  Satis¬ 
fying  global  boundary  conditions  affects  the  steady  pressure 
level,  about  which  the  dynamic  motion  of  the  damper  causes 
a  perturbation.  Satisfaction  of  local  boundary  conditions 
Influence  the  dynamic  pressure  response  by  Influencing  the 
local  dynamic  continuity. 

Frequency  Is  accounted  for  in  a  similar  manner  to  that 
already  described  for  finite  velocity.  Indeed,  whether 
one  is  concerned  about  the  effect  of  finite  velocity  or 
finite  frequency  is  dependent  upon  whether  the  general 
impedance  analysis  is  being  performed  or  the  more  specialized 
circular  orbit  analysis.  In  the  case  of  the  circular  orbit 
analysis,  the  orbit  frequency  sets  the  velocity.  For  use 
in  rotordynamics  calculation,  the  frequency  is  the  most 
convenient  parameter  to  talk  about.  For  fluid  film  compu¬ 
tation,  the  implicit  finite  velocity  is  accounted  for  and, 
in  combination  with  supply  conditions,  dictates  the  loca¬ 
tion  of  the  zero  pressure  boundaries,  and  thereby,  the 
extent  of  the  cavitated  region. 

The  direct  result  of  the  Reynolds  equation  solution  are 
pressures.  These  are  integrated  numerically  to  give 
forces.  The  relationship  of  force  to  velocity  or  dis¬ 
placement  gives  coefficients  as  needed. 

Flow  is  calculated  by  integrating  around  the  boundary  of 
interest.  The  local  flow  is  determined  from  the  computed 
pressure  gradient,  film  thickness,  and  viscosity. 


The  following  is  a  description  of  the  damper  film  analysis,  along  with  des¬ 
criptions  of  the  solution  scheme,  the  mobility  method  application  and  jump 
analysis. 
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2.1  Reynolds'  Equation  Without  Inertia  Effects 

The  squeeze  film  damper  Is  a  cylindrical  bearing  which  supports  a  non¬ 
rotating  bushing  or  journal.  Figure  2-1  shows  the  fixed  coordinate  system 
and  quantities  used  to  describe  the  damper  geometry. 


x 

Figure  2-1  Squeeze  Film  Damper  Geometry 

The  bearing  dimensions  and  lubricant  properties  are: 

D  *  bearing  diameter,  (m) 

R  *  1/2D  ■  bearing  radius,  (m) 

L  -  bearing  length,  (m) 

C  «  radial  clearance,  (m) 

2 

U  ■  lubricant  viscosity, (N-sec/m  ) 

3 

p  ■  lubricant  density, (kg/m  ) 
x,  y  ■  bushing  center  coordinates,  (m) 

9  *  circumferential  coordinate  angle,  measured  from  negative  x-axls  (radians) 
z  ■  axial  coordinate,  (m) 

h  ■  C  +  xcos0  +  ysin9  ■  local  film  thickness,  (in)  (2-1) 
t  -  time,  (sec) 

*  flow  from  k'th  feeding  hole,  (m  /sec) 
f  m  reference  frequency,  (Hz) 

u  ■  2xf  ■  reference  angular  frequency,  (radians/sec) 

2 

p  •  film  pressure,  (N/m  ) 
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By  ignoring  inertia  effects,  the  pressure  p  in  the  lubricant  film  derives 
from  Reynold's  equation  in  the  fora: 


3  3 

_L_  (  h  9P\  +  _JL/_JL. 
R39  '12y  R39;  3zv12y 


1£) 

3z; 


lii  I 

3t  "  R 


l  q^  6(9  -  9^)  d(z  -  z^) 


(2-2) 


where  6(9  -  9^)  and  S(z  -  z^)  are  unit  delta-functions,  and  9^  and  z^  give 
the  coordinates  of  the  k'th  feeding  hole.  Furthermore: 


3h 

3t 


dT  C080  +  at  sin9 


(2-3) 


With  z  measured  from  the  centerplane  of  the  bearing  the  boundary  conditions 
for  eq.  2-2  are: 

at  z  ■  1/2L:  '  2*  3 

with  sealed  end  groove:  I  (-  f^i/2L  1140  "  CG  P^^L)  (2-4) 


with  end  seal  (no  groove) : 


JL.  !£■) 

^  12y  3z'l/2L 


-  p(l/2L)  (2-5) 


with  no  end  sealing:  p(l/2L)  •  0  (2-6) 

where  Cg  and  <3^  are  flow  coefficients  for  the  seals.*  When  the  bearing  is 
not  symmetric  with  respect  to  the  centerplane,  these  boundary  conditions 
also  apply  at  z  ■  -  -|l.  When  there  is  symmetry  around  the  centerplane,  only 
half  the  film  needs  to  be  considered,  and  the  boundary  conditions  at  z  ■  0 
become : 


at  z  -  0: 


with  circumferential  feeding  groove:  p(0)  ■  p 


(2-7) 


without  circumferential  feeding  groove:  (ff}0  "  0  <2'8> 


where  pg  is  the  supply  pressure  in  the  groove. 


*  Flow  coefficients  for  a  resistance  relate  flow  through  the  resistance  to 
pressure  drop  across  the  resistance  units  are  (m^/N’sec)  for 
(m^/N*sec)  for  Cl*  ^ 
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The  pressure,  pk»  at  the  k'th  feeding  hole  is  related  to  the  feeder  hole 
supply  pressure,  pgk>  and  the  flow,  q^,  by: 

at  a  -  ak,  .8  -  9k:  Pk  -  P,k  -  ^  9k  (2-9) 

where  represents  the  resistance  of  the  feeder  hole. 


In  order  to  generalize  the  results  the  equations  are 
form.  The  variables  are  made  dimensionless  by: 


t  -  a)  t 

p  -  p/12um(^)2 
5k-  qk/uCR2 
x,  y  -  x/C,  y/C 

h  -  h/C  -  1  +  x  cos©  +  y  sin© 


solved  in  dimensionless 


(2-10) 

(2-11) 

(2-12) 

(2-13) 

(2-14) 

(2-15) 


Thereby  eq.  (2-2)  is  transformed  into: 

•35  (h3  +  -3“(h3  J^)  •  x  cos©  +  y  sin©  -  S  qfc  6  (0 

where 

-  -  3x  _i  dx  _i  dj. 

x*  y  "  TF*  3t  “  Cm  dt’  Cu  dt 


©k)  «(C  -  Ck)  (2-16) 


(2-17) 


a>  (radians/sec)  is  a  representative  angular  frequency.  Its  precise  value  is 
not  important  to  the  dimensionless  form  of  the  equation  as  it  only  serves  as  a 
scaling  factor,  but  typically  it  should  be  selected  such  that  x,  y  are  of  the 
order  of  1. 

at  C  -  £  (*  -  -jL) : 

with  sealed  end  groove: 


2w 

/' 


(-  h3 


12.) 

3?;L/D 


d©  ■  C  p(C  ■  L/D) 


(2-18) 
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(2-19) 


with  end  seal  (no  groove):  (-h  "  p) 


with  no  end  sealing:  p(5  ■  L/D)  ■  0 


(2-20) 


The  dimensionless  seal  flow  coefficients  are  defined  as: 
12u 


C_  «  •=#=■  C-  (sealed  end  groove) 

G  c3  G 

C  .  c  (end  seal,  no  groove) 

L  C3 * S 


(2-21) 

(2-22) 


As  an  example,  the  seal  may  be  represented  as  a  ring  of  width  W  with  a 
leakage  gap  of  height  S,  as  shown  in  Figure  2-2. 


wvwww 


£ 


.\\\\ 


w 


ambient  pressure  •  0 


r 


Figure  2-2  Representation  of  Sealing  Ring  in  Groove 


The  flow  through  the  seal  gap  is: 

3 

for  sealed  end  groove:  2xR  P »  Cg  p(l/2L)  (see  Eq.  (2-4)) 

3 

for  end  seal,  no  groove:  p(l/2L)  (see  Eq.  (2-5)) 

With  this  seal  model  the  dimensionless  seal  coefficients  become: 

CG  -  2*  (£)  (|)3  (2-23) 

CL  -  (§)  (f)3  <2-24) 

When  there  is  symmetry  around  the  centerplane,  the  boundary  conditions  at 

S  ■  0  are: 


11 


at  ?  ■  0: 


with  circumferential  feeding  groove:  p(0)  «  p 


without  circumferential  feeding  groove:  Gr?-)  -  0 

o 


(2-25) 

(2-26) 


where  pg  is  -the  dimensionless  supply  pressure: 


ps  "  12uu»(|)2 


(2-27) 


At  the  feeder  holes,  Eq.  (2-9)  applies: 


“  5  ■  V  9  •  V  *k  ■  -  \  \ 


(2-28) 


where  is  the  dimensionless  feeder  hole  resistance: 
C3 

°  k  "  I2y  “k 


(2-29) 


As  an  example,  the  feeder  hole  may  be  represented  as  a  capillary  for  which  the 
flow  is  given  by: 

,4 


2:rd 

Vi 


<pak  -  pk) 


(2-30) 


where  d  is  the  tube  diameter  and  £  is  the  length.  Hence,  the  resistance  is: 
ul 

S' 77 


and  the  dimensionless  resistance  becomes: 


-  _  1  ,1.  ,C.3 

ak  *  247  (d>  (d} 


With  the  given  boundary  conditions  Eq.  (2-16)  may  be  solved  numerically  as 
discussed  later. 

By  integrating  the  film  pressures  the  damper  film  forces  are  computed  as 
(see  Fig.  2-1): 


12 


2ir 


I 


COSGI 


sinOl 


R  d9  dz 


or. 


In  dimensionless 

r - SL-7 

*  u«ta(|)2 

Fy 

F  -  .R.  2 

y  WuLR(-^) 


form: 

6 

►  "  L/D 


(2-31) 


(2-32) 


The  flow  out  of  the  film  at  the  end  is: 


at  z 


(-  _h^  3£)  RdQ 

12p  3z  -k 


(2-33) 


which  in  dimensionless  form  becomes: 


at  z 


L 

D' 


/  <_i'3  h’l/d  d0 


O 


(2-34) 


'  L 

An  analogous  expression  holds  for  the  end  at  f  -  and  for  the  flow  in  from  a 
circumferential  groove.  The  inflow  from  any  feeder  holes  is  the  sum  of  the 
calculated  values. 


2.2  Inertia  Effects  in  Reynold’s  Equation 

Under  the  conventional  assumption  of  lubrication  theory  where  the  lubricant 
film  is  assumed  thin  (C/R<<1),  the  Navier-Stokes  equations  reduce  to: 


P<!r  +  ui^+v!7+w!7> 


0(|2i+ui2_.  +  vis  +  w.i2)  .-i£ 
V  u  R3Q  3.y  3z'  3z 


(2-35) 
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where  y  is  the  coordinate  across  the  film,  u,  v  and  w  are  the  fluid  velocity 
components  in  the  9,  y  and  z  directions,  respectively,  and  0  is  the  fluid 
mass  density.  The  terms  on  the  left  hand  side  of  the  equations  represent  the 
inertia  forces,  and  they  are  normally  ignored  as  being  small  in  most  cases. 

To  assess  their  magnitude,  the  equations  are  made  dimensionless  by: 

-  u  -  w 

u  ■  —  w  -  — 

Rio  Rid 


v 

Cio 


-  L 

C 


while  the  remaining  variables  are  made  dimensionless  in  accordance  with  eqs. 
(2-10),  (2-11),  (2-12)  and  (2-15).  The  dimensionless  equations  become: 


'"■Tt 


(2-36) 


u 


3w 

39 


-  3W  -  3w. 
+  v  3^  +  »  T? 


where 


and  Re  is  the  Reynold's  number: 
r 

e  u 


(2-37) 


(2-38) 


Thus,  only  when  X  is  of  the  order  of  unity  are  the  inertia  terms  important. 

To  include  them,  a  perturbation  solution  is  employed,  settins: 

p  «  pQ  +  X  px  +  — 

(2-39) 

u  ■  u  +  X  u,  +  — 

•  o  1 

and  similarly  for  v  and  w.  By  substituting  into  Eqs.  (2-36)  and  collecting  terms 
of  like  power  in  X,  the  following  equations  result: 


(2-40) 
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(2-41) 


3  u.  3p  3u  3u  3u  3u 

3n^  "  12  ~~ 3®  +  TF  +  uo  ~39  +  Vo  In’  +  Wo  Tc 


(2-42) 


}  W1  .  3pl  .  3wo  ^  -  3wo  .  -  3wo  A  -  3wo 
TP  12  Tc  +  It  u0  Te  +  vo  TrT  wo  Ts 


(2-43) 


(the  equations  for  higher  order  terras  are  Ignored) , 


As  the  pressure  is  independent  of  jr,  p  and  p,  do  not  depend  on  n,  and 

o  i 

eqs.  (2-40)  and  (2-41)  may  be  integrated  twice  to  give: 


Uo  ’  ‘  6  T©  n(h  “  n) 


(2-44) 


Wo  "  “  6  Tc  n(S  “  n) 

which  satisfy  the  boundary  conditions  that  u  and  w  are  zero  at  the  two  film 
surfaces  (y  »  0  and  y  -  h) .  The  continuity  equation  is 


3u  .  3v  ,  3w  . 
+  3^+  3?"  0 


(2-45) 


-  -  2  -  2- 
3u  3w  3  p  _  3p  r  3  pa 

a  -  it  - 6  — ? 11  (h  -  +  6 1?  n " + 6  tt  -  *» 

39  ®  4 


and  upon  integration: 

2-  2- 
_  O  3  P_  j  it  -y  3  P_ 

vQ  •  n  (3h  -  2n)  y^r  +  3n  §-§■  +.n(3h  -  2n)  - y 

34 

(boundary  condition:  t  ■  o  it  n  ■  o). 


(2-46) 
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At  n  ■  h,  v  equals  the  squeeze  velocity  «-=■  such  that  Eq.  (2-46)  yields: 

O  3T 


(2-47) 


which  is  Reynold's  equation  without  inertia  effects  as  previously  given  in 

Eq.  (2-16).  By  substituting  into  Eq.  (2-46),  v  may  also  be  expressed  as: 

o 


n2(3h  -  2n)  dh  6n2(h  -  n)  3h  3po 

Vo  *  g3  dt  g  39  39 


39  39 


(2-48) 


With  u  ,  w  and  v  determined,  Eqs.  (2-41)  and  (2-42)  can  be  integrated  twice  to 
o  o  o 

yield: 

2-  - 

-  -  ^  P1  1  -3  “3  4  ^  Po  -2  2  3po 

-  -6n(h  -  n)  +  Y(nh  -  2hnJ  +  O  +  nOT  -  n  )  ■£ 


■//<=, 


3u  3u  3u 

:0lf+’o?r+i;o  “5?)dndn  (2-49) 


—  2-  _ 
wx  -  -6n(h  -  n)  -|i.  +  I(nh3  -  2hn3  +  n4)  -gjf  +  n(h2  -  n2) 


ff  3w  ^  3w  3w 

'JJ  (Zo  “4  +  *o  "4  +  *o  1?)dndn  '  ~50) 


By  restricting  the  analysis  to  the  terms  which  Involve  the  squeeze  acceleration 
2- 

— ■  x  cos9  +  y  sin9,  only  the  two  first  terms  on  the  right  hand  side  of  the 
dt 

equations  are  retained,  after  which  u^  and  are  substituted  into  the 
continuity  equation: 


such  Chat: 


3v. 


_1 

3n 


a2;, 


.3- 


6n(h  -  n)  |  +  6n  f|  -1  -  \  (nh3  -  2hn3  +  -i4)  -~^2- 

3Q2  393Q-  2  3923T 

i  *)  sc  ^  p  d  p. 

-i(3nh2-2n3)5|j^f+6n(5-,)-J 

-  i  <„h3  -  2hn3  +  „4)  i^2- 

30x 


This  equation  may  be  Integrated  with  respect  to  n,  applying  the  boundary 
conditions  of  v  »  0  at  n  -  0  and  n  -  h.  By  making  use  of  Eq.  (2-47)  and 
retaining  only  acceleration  terms,  the  final  result  is: 


<h  ^e>  a5)  10  3^2  5  39  3739 


(2-51) 


Trom  Eqs.  (2-47)  and  the  boundary  conditions  it  is  seen  that  pq  may  be 
written  as  the  sum  of  3  terms: 


po  ■  p«  *  +  Of  1  +  pb  P, 


(2-52) 


Similarly  p^  may  be  written  as: 


Pi  -  Pk-  x  +  p?  y 

By  substitution  into  Eq.  (2-51)  two  equations  are  obtained: 


-3  3px 
J _ *\ 


-3  3Px, 


f?<h  +  rt> 


3  _■*  3P»  a  .1  3Pu 

39<h  +  T?<h  7?*- 


~  coafl  +  11.  lE  9P* 

10  C0*e  5  39  39 


£  .me  +  £  !£  !li 
To  *in®  +  ~  Is  "Tq 


(2-53) 


(2-54) 
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2.3  Solution  of  Reynold's  Equation 


E 

I 


•1 

j 

i 

1 


i 

< 


The  solution  of  Eqs.  (2-47)  and  (2-16)  is  obtained  in  the  form: 


_  „i  i  i  _ 

P  *  p  .x  +  p  .  y  +  Sp  q_  +  p  .  p 
ro  r  x  ■  y  o  mm  brs 


where  the  variables  are  solutions  of: 


(i«<53  y + 1?)} 


(2-55) 


p  • 

X 

cos© 

(a) 

—  T 

P  • 

sin© 

(b) 

7 

v  -  i 

1  1 

P  - 

-6(9  -  0)  5(5  -  C) 

(c) 

m 

y 

m  m 

P  V 

0 

(d) 

L  “J 

J 

\  (2-56) 


The  equation  is  solved  as  a  finite  difference  equation  as  discussed  later. 


The  boundary  conditions  are  as  follows:  if  there  is  end  seal,  but  not  in  a 
groove,  Eq.  (2-19)  applies  such  that: 


at  5 


L 

D* 


CL(P  x^L/D 


(2-57) 


»  i  *  i 

and  similarly  for  p  p  and  p  ,  .  The  solution  for  p  .  ,  however,  is  only 

y  tn  t>  0 

needed  when  there  is  a  circumferential  feeding  groove.  If  the  bearing  is  not 
synanetrlc  around  the  centerplane,  Eq.  (2-57)  also  applies  at  ?  *  L/D.  If  there 
is  symmetry,  the  boundary  conditions  at  the  centerplane  are: 


at  C  ■  0: 


with  circumferential  feeding  groove: 


p  .  »  p  .  -  p  ■  0,  p  .  ■  1 

x  y  m  rb 


no  circumferential  feeding  groove: 


I  f 

dp  .  dp 

d?  "  dT 


0 


(2-58) 


t 

' 
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If,  instead,  the  end  seal  is  in  a  groove,  the  boundary  condition  at  the  end 
is  given  by  Eq.  (2-18).  The  groove  pressure  may  be  expressed  as: 


p0(c  “  b  “  p*+)* +  H+)? +  z  +  Pb+)ps 


(2-59) 


If  the  bearing  is  not  symmetric  around  the  centerplane,  the  groove  pressure 
at  the  other  end  is:  PQ(C  ■  -*^)  *  Pq  ^  which  may  be  expressed  analogously. 

The  pressure  field  set  up  by  the  groove  pressure  and/or  the  pressure  in  the 
circumferential  feeding  groove  is: 


(1  -  -j4j))  Ps  +  (with  circ.  feed  groove,  symmetric 

'  bearing) 

(2-60) 

p^+^  (without  circ.  feed  groove,  symmetric 

'  bearing) 


yd  -  P^  ^  +  -j(l  +  ^7^)  p^+)  (not  symmetric  bearing) 


which  replaces  the  term:  p  .  p  in  Eq.  (2-55)  when  the  end  sea.  are  in  a  groo 

0  s 

By  substitution  into  Eq.  (2-18),  the  boundary  :-;«uiitiotitt  at  q  «  l/D  becomes: 


-3  P  x  - 

(-h  -jf)  dQ  +  y 


tv 


-3  3p  y 

»  -H1'  d9  + 


4m  " 

Z  \  / 
m  4, 


3P  » 

(-b  d9 


L.  » 

I 


53  d9  •  C’ 


(2-61) 


When  the  film  is  fed  from  a  circumferential  groove,  and  not  from  feeder  holes 
such  that  the  q_  term  does  not  apply,  the  equation  can  be  solved  directly  for 
p^  to  get  the  pressure  parameters  in  Eq.  (2-64). 
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On  the  other  hand,  when  there  are  feeder  holes  and  no  circumferential  feeder 
groove,  Eq.  (2-61)  must  be  solved  simultaneously  with  the  feeder  hole  flow 
equations.  They  are  obtained  from  Eq.  (2-28)  with  substitution  from  Eqs. 
(2-55)  and  (2-60): 

I  _  »  .  • 

at  C  “  V  9  "  V  ®1c  +  J  p  a  in  +  Pb  ~  Pyy  +  ^sk  (2~62) 

These  equations  can  be  solved  together  with  Eq.  (2-61)  to  give  q,  and  p^ 

-(-)  k  ° 

(and  p  )  in  the  form: 
o 


A  -  - 


"k  *  ’kx  x  +  V  y  +  "VB  P» 


-(+)  -(+)  *  ,  -(+)  *  ,  -(+)  - 
P„  “  p  .  x  +  p y  +  p__  p 
o  ox 


p  .  y  +  p  p 
oy  os  *s 


(2-63) 

(2-64) 


By  substituting  these  equations  into  Eqs.  (2-55)  and  (2-60),  the  pressure  distri¬ 
bution,  p  ,  can  be  expressed  as  given  in  Eq.  (2-52).  By  integrating  p  in 
o  o 

accordance  with  Eq.  (2-32),  the  dimensionless  film  forces  are  obtained  in 
the  form: 


F  *  -  B  x  -  B _ y  +  F 

x  xx  xf  xs 


F  •  -  B  x  -  B _ y  +  F 

y  yx  rr  ys 


(2-65) 


where: 


"HhO 

-T  « 


d0dc 


(2-66) 


£  2* 


6 

l7d 


If 


'cose" 


sin©. 


dOdc 


(2-67) 
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and  B 

yy 

B 

XX 


UL(|)3 


and  similarly  for  Che 


are  che  dimensionless  damping  coefficients  in 


ocher  coefficients. 


Once  p.  and  p.  are  determined,  Eq.  (2-54)  can  be  solved  for  p-  and 
x  y  x 

can  be  integrated  in  accordance  with  Eq.  (2-32)  where  p  is  replaced 

Xp  »  MP"X  +  p..y)  (see  Eqs.  (2-39)  and  (2-53)): 
j.  x  y 


(2-68) 

the  form: 

(2-69) 

which 

y 

by 

(2-70) 

(2-71) 


(2-72) 
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C  ,  C  and  C  are  the  dimensionless  inertia  coefficients  in  the  form: 

xy’  y  x  yy 


(2-73) 


and  similarly  for  the  other  coefficients.  The  actual  coefficients  have  the 
unit  of  mass. 


2.4  Numerical  Solution 

Reynolds'  equation  is  solved  as  a  finite  difference  equation.  A  finite 
difference  grid  shown  in  Figure  2-3,  is  introduced  with  a  j-axis  in  the 
circumferential  direction  (9-direction)  and  an  i-axis  in  the  axial  direction 
(z-direction) : 


j-0  j-1  j-2 -  j«m 


Figure  2-3  Finite  Difference  Grid  for  Damper  Film 


When  the  bearing  end  is  a  sealed  groove,  or  there  is  no  seal,  the  end  is  at 
i  ■  0.  When  the  seal  has  no  groove,  the  end  is  at  i  •  1.  For  a  symmetric 
bearing  where  only  one  half  of  the  film  needs  to  be  considered ,  the  center- 
plane  is  at  i  -  n.  For  a  full  bearing,  i  -  n  gives  the  other  end  of  the  film. 


The  finite  difference  increments  are: 


A/i  2ir 

49  ■  — 

IS 


(2-74) 
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wh«re  a'  is  adjusted  according  to  whether  the  film  is  half  or  full,  and  to  the 
sealing  conditions. 

In  finite  difference  form,  the  dimensionless  Reynold's  equation  (Eqs.  (2-56) 
and  (2-54))  may  be  written  as: 


bij  pi,j-i +  6i.j«  +  *'«  h-1,3  *  *u  hi  *  *"«  h+i,i  -  f 

(2-75; 

■trtiere  p  represents  any  of  the  six  pressure  variables  in  Eqs.  (2-56)  and  (2-54), 
and  where: 


b  -  — i— 

«  de2  lb^1/2 


b  -  -L-h3 
°i,j+l  lQ2  ni,j+l/2 


(2-76) 


.J-h3 

ij  4^  ni-l/2, j 


1:3 

ij  li?  i+l/2,j 


(2-77) 


“u  *  -  (*ij  +  *«  +  b«  +  bi,j+i> 


(2-78) 


f^  *  right  hand  side  of  Eqs.  (2-56)  and  (2-54),  evaluated  at 


9  ■  »j  .  '  -  5i 


For  the  special  case  of  Eq.  (2-56c) ,  f^  is  given  bv 


ij 


0  when  9.  5*  0  or;,  i*  ; 

j  m  i  m 


when  -  9n,  ^  ^ 


(2-79) 


When  the  bearing  end  has  a  sealed  groove,  or  there  is  no  seal,  then: 


poj  *° 


(2-80) 


When  the  bearing  end  has  a  seal  which  is  not  in  a  groove,  Eq.  (2-19)  yields: 

f  Pij  +  Z  P^4  -f  PoW  ■  C,  P,,  (2-81) 


E3 


4^  <-  p2j  -}p3j)  ■  tj.  hi 


4? 
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where  the  gradient:  -  ^  is  derived  on  the  basis  of  a  second  order  polynomial. 

dC 

By  comparing  with  Eq.  (2-75)  it  is  found  that  in  this  case: 


i  -  1  : 


bij  "  \  j+1  “  fiJ 


0 


K3 

+ 1  ii* 

2i? 


> 


2^- 
2AC  p2j 


0 


(2-82) 


For  a  full  length  film,  the  boundary  conditions  at  the  opposite  end  (i  *  n) 
are  analogous  to  Eqs.  (2-80)  and  (2-82).  For  a  symmetric  fi*m,  i  ■  n  gives 
the  centerplane  where  the  boundary  conditions  are  given  by  Eqs.  (2-25)  and 
(2-26)  : 

at  i  -  n:  with  circumferential  feeding  groove:  ■  0  (2-83) 

without  circumferential  feeding  groove:  p  , ,  .  »  p  .  . 

n+l,J  n-i,J 


The  pressures  along  any  j-grid  line  may  be  arranged  in  a  column  vector: 


Ti 


whereby  Eq.  (2-75)  can  he  written  as  a  matrix  equation: 


(2-84) 


PJ  +  Pj-1  +  ~)+l  ?j+l  "  Fj 


where  F^  is  the  column  vector: 

rt  ' 
lj 


F*  -  <£ij> 


(2-85) 


(2-86 
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Bj  and  B^+^  are  diagonal  matrices: 


%  -  dl**  {bli>  1,4.1  '  dl**  {bl.J+l) 


(2-87 


and  Aj  is  a  tri-diagonal  matrix: 


* 
r 


\  <9* 

\  \ 

<9  A 


\ 


\ 


im  (  y  ^ 

'  V  \  ♦* 

\  o  A 

\ 

tr 


(2-88 


It  should  be  noted  that  the  boundary  conditions  of  eqs.  (2-80),  (2-82)  and 
(2-83)  are  built  in  A ^ and  B^+^.  Eq.  (2-85)  is  solved  in  the  form: 

(2-89 


P.  1  "  1  +  <L  ,  P  +  E.  . 

j-1  “j-1  j  -j-1  m  j-1 


where  and  are  matrices  while  is  a  column  vector.  By  substitution 
into  Eq .  (85)  it  is  found  that: 


D  ■  -  H  B 

^  Vi 


G,  -  -  H,  B  G,  , 

-j  -j  -i  -j-i 


(2-90 


WWi-i’ 


where 


&  -  (A,  +  B. 


-1 


r)  xsj  -=J  =j-l 


(2-91 


As  the  pressure  at  0  ■  0  has  to  be  equal  to  the  pressure  at  9  -  2ir,  circum¬ 
ferential  continuity  requires  that: 


Hence,  it  Is  seen  from  Eq.  (2-89)  that: 


^  -  0  Eq  -  0  (^  -  I  (2-93) 

where  _I  is  the  identity  matrix.  Hence,  Eqs.(2-90)  and  (2-91)  may  be  used  for 
j  going  from  1  to  m.  For  j  »  o+l,  eq.  (2-89)  yields: 


P  -  D  P  . ,  +  G  P  +  E 
m  — m  ofl  — m  — m  m 


With 


P^^  •  this  equation  may  be  solved: 


-1 


Pm  "  P„  ■  (I  “  6  )  *  (D  P  +  E  ) 
m  o  -m  — m  l  m 


(2-94) 


With  P^  still  unknown,  set: 


Wi + 


(2-95) 


where  S.  is  a  matrix  and  T  is  a  column  vector.  It  is  seen  from  Eq.  (2-94) 
that:^  j 


(I  -  G  r1  D 


T  -  (I  -  G  r1  E 

uX  u2  m 


(2-96) 


By  substituting  Eq.  (2-95)  into  Eq.  (2-89)  it  is  found  that: 


•=j-l  ~j-i  -^j-l  -a 


T.  i  *  D.  ,  T,  +  G.  .  T  +  E,  . 
j-1  -j-1  j  -j-1  m  j-1 


(2-97) 


With  S  and  T  given  by  Eq.  (2-96),  Eq.  (2-97)  may  be  used  for  j  going  from 
— m  tn 

J  -  m  to  j  -  2  after  which  Eq.  (2-95)  yields: 


(2-98) 
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(2-99) 


This  equation  may  be  solved  for  P^: 


after  which  Eq.  (2-95)  can  be  used  to  compute  P^  for  j  going  from  j-2  to 
j-m  whereby  the  complete  pressure  distribution  is  determined. 


The  advantage  of  this  method  is  that. the  most  time  consuming  part  of  the 
calculation,  namely  the  calculation  of  (Eq.  (2-81))  for  j  going  from  j*l 
to  j*m,  needs  only  to  be  done  once  in  any  given  case.  Thus,  Eqs.  (2-90), 
(2-91),  (2-96)  and  (2-97)  may  be  employed  to  compute  and  store  and 

(and  also  (I  -  G  )_1  and  (I  -  SL)”1,  see  Eqs.  (2-96)  and  (2-99)).  With  B, 

—  ""TO  “1  —J 

stored  during  the  calculations,  Eqs.  (2-90),  (2-96),  (2-97)  and  (2-99)  and, 
finally,  Eq.  (2-98),  may  be  used  to  obtain  the  solution  for  any  of  the  six 
pressure  variables  in  Eqs.  (2-56)  and  (2-54)  by  simple  matrix  multiplication. 
The  stored  matrices  represent  the  operator  on  the  left  hand  side  of  Eqs. 
(2-56)  and  (2-54)  together  with  the  associated  boundary  conditions  which 
means  that  they  apply  at  a  given  journal  center  position,  a  specified  seal 
leakage  coefficient  for  an  end  seal  without  a  groove,  and  whether  the 
bearing  is  symmetric  or  of  full  length.  To  make  calculations  for  any  number 
and  combinations  of  feeder  holes,  or  sealed  end  grooves  with  difference 
leakage  coefficient,  requires  little  computation  time.  The  boundary  condi¬ 
tions  to  satisfy  are  those  of  Eqs.  (2-67)  and/or  Eq.  (2-62). 


The  calculations  end  up  with  five  pressure  distributions,  namely  p.,  p.  and 

__  _  x  Y 

(see  Eq.  (2-52)),  and  p^  and  p^.  (see  Eq.  (2-53)).  They  are  integrated 
according  to  Eqs.  (2-66),  (2-67),  (2-68),  (2-71)  and  (2-72)  to  yield  the 


four  damping  coefficients  and  the  four 
are  evaluated  numerically  as: 


inertia  coefficients.  The  integrals 
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t  a 


(The  result  Is  multiplied  by  two  when  only  half  of  the  film  has  been  computed. 
Also,  any  point  at  the  end  or  at  the  centerplane  only  contributes  by  half  its 
pressure  value) . 

The  integration  is  performed  for  given  values  of  the  squeeze  velocity  (x  and  y) 
and  of  the  supply  pressure,  pg.  -Hence,  the  pressure,  pQ^  can  be  computed  at  any 

selected  point  from  Eq.  (2-52).  If  the  thus  calculated  value  is  less  than  a 
specified  cavitation  pressure,  the  pressure  at  the  point  is  set  equal  to  the 
cavitation  pressure,  and  the  point  does  not  contribute  to  the  damping  coeffi¬ 
cients  or  the  inertia  coefficients.  This  means  that  only  the  pressures  in  the 
non-cavitating  domain  of  the  film  are  included  in  the  integration,  and  as  the 
domain  changes  with  the  magnitude  and  direction  of  the  squeeze  velocity,  so  do 
the  film  forces  and  coefficients  change.  This  also  means,  that  the  calculated 
damping  coefficients  are  not  exactly  equal  to  the  partial  derivatives  of  the 
film  forces. 


2.5  Comments  on  the  Mobility  Method 

The  external  forces  acting  on  the  damper  journal  are  W  and  W  ,  shown  in 

x  y 

Figure  2-4,  which  are  functions  of  time. 


Figure  2-4  Forces  Acting  on  the  Journal 


The  forces  from  the  film  are  F  and  F  which  are  function  of  the  journal 

x  y 

center  coordinates,  x  and  y,  and  the  corresponding  velocities  x  and  y.  Let 
the  journal  at  time  t  be  at  the  position  x(t) ,  y(t) .  A  force  balance  requires 
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W,  +  p-  -o 


Wy  +  Fy  -0 


(2-101) 


These  equations  are  satisfied  by  determining  those  values  a*t  x  and  y  which  will 
make  the  expressions  zero.  Thus,  with  guessed  values  at  x  and  y  an  expansion 
yields: 

■  / 3Fx \  Ax  /3FX|  Ay 

Wx  +  Fxo  +  |3i|0  +|3^"j0  -0  (2-102) 


(2-102) 


and  similarly  for  the  y  direction.  Fxo,  and  the  derivatives  are  evaluated 

for  the  given  values  at  x(t),  y(t),  x  and  y.  Set: 

/3F„l  /  3F  \ 


B_  -  - 


O 


_ X 

3^  ,  o 


yx  \  3j  fo 

such  that  eq.  (2-102)  nay  be  written: 


3  F 

- - z 

t*n« 


(2-103) 


W  +  F  > 
X  xo 


W  +  F 
•  7  71 


(2-104) 


This  equation  is  readily  solved  for  Ax  and  Ay  .  Because  B  ,  B  ,  B  and  B 

xx  xy  yx  yy 

in  themselves  are  functions  of  velocity,  due  to  the  change  in  cavitation  boundary, 
it  is  necessary  to  iterate  a  number  of  times  until  the  right  hand  side  becomes 
sufficiently  small  and  thereby  satisfies  eq.  (2-101). 

•  • 

With  x  and  y  determined  at  time  t,  the  next  journal  center  position  is: 


x(t  +At)  -  x(t)  +  Afx(t)  y(t  +At)  -  y(t)  +&t*y(t)  (2-105) 

after  which  the  procedure  is  repeated. 

Fx,  Fy  and  the  four  damping  coefficients  are  computed  by  interpolation  from 
stored  pressure  field  data.  These  data  assume  that  the  bearing  geometry  is 
axisymmetric  such  that  the  only  parameter  is  the  eccentricity  e  as 
illustrated  in  Figure  2-5. 
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Brg.  Ctr.-^ 
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Figure  2-5  Forces  in  Sri  Coordinate  System 


Hence,  the  stored  data  refers  to  the  £n  coordinate  system,  and  can  be  used 
to  compute  F^,  and  the  4  associated  damping  coefficients  ,  B^ 

and  B^n  .  To  transform  back  into  the  fixed  x-y-system  it  is  seen  that: 


e  * 


/7T? 


COSO  “  ^ 


sin  a 


(2-106) 


and: 


F_ 


Fr  cosa  -  F  sin® 

£  n 


F  ■  F_  sina  +  F„  cosa 
y  4  n 


5  »  x  cos  a  +  y  sina 

while  the  coefficients  transform  as: 
B 
B. 

B 
B 


n  -  -  x  sina  +  y  cos  a 


(2-107) 
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(2-108) 


30 


2.6  Static  Jump 


2.6.1  Analysis  Using  Short  Bearing  Approximation,  Rigid  Rotor 

Static  jump  is  a  condition  wherein,  for  a  given  set  of  conditions,  a  damper 
journal  may  have  more  than  one  stable  orbit.  In  fact  as  many  as  three 
possible  orbits  may  exist:  a  "whirl  orbit"  with  large  eccentricity,  an 
"inverted  orbit"  with  small  eccentricity,  or  an  intermediate  orbit  which 
is  not  stable  and  degenerates  into  one  or  the  other  of  the  first  two.  The 
shift  from  one  steady-state  orbit  to  the  other  is  discontinuous,  hence  the 
use  of  the  term  "static  jump"  to  describe  the  change  between  states.  White, 
in  Reference  (2)  performed  an  analysis  of  static  jump  in  which  he  used 
the  "short  bearing"  -approximation  for  the  squeeze  film  and  certain  force 
relationships  to  derive  charts  or  jump  "maps".  The  procedure  described 
below  is  similar  to  White's  but  extends  the  analysis  to  finite  length 
dampers  with  parallel  springs  and  rotor  flexibility. 

The  first  problem  to  be  addressed  is  to  develop  the  conditions  for  equili¬ 
brium  (i.e.,  a  stable  orbit)  for  a  whirling  damper  executing  a  circular 
orbit  to  support  a  rigid  rotor  whose  mass  center  is  eccentric  from  the 
geometric  center  of  the  damper  by  a  distance  etc,  where  c  is  the  damper  radial 
clearance  and  a  is  a  nondlmensional  unbalance  parameter.  Figure  2-6  shows  the 
various  quantities  involved  in  the  analysis  and  in  particular  3  points:  the 
center  of  the  damper  clearance  circle,  the  damper  geometric  cer.cer  and  the 
mass  center  of  the  rotor.  The  distance  between  the  mass  center  and  the 
geometric  center  of  the  damper  is  defined  by  the  unbalance  parameter  a  and 
remains  fixed  for  a  particular  equilibrium  analysis.  If  the  angular  velocity 
of  the  rotor  is  w,  a  synchronous  whirl  orbit,  also  with  frequency  w,  will 
be  developed. 

As  shown  by  Figure  2-6  two  distances  define  the  orbit  which  is  developed. 

The  first  is  the  distance  from  the  center  of  the  clearance  circle  to  the 
geometric  center  of  the  damper.  This  distance  is  expressed  as  the  product 
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Figure  2-6  Equilibrium  Geometry  for  Whirling  Damper  Supporting  a  Rigid  Rotor 
with  Mass  Center  Eccentric  from  Geometric  Center  by  Distance  ac 
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of  e  and  c,  where  e  is  a  uond linens ional  orbit  circle  radius  and  c  is  the 
radial  clearance  of  the  damper.  The  second  is  the  distance  from  the  center 
of  the  clearance  circle  to  the  center  of  mass.  This  distance  is  expressed 
as  the  product  of  3  and  c,  where  6  is  a  nondimens ional  mass  center  orbit 
radius. 

To  quantify  the  severity  of  whirling  vibration  exhibited  by  this  unbalanced 
rotor  in  a  particular  squeeze  film  damper,  the  problem,  mathematically,  is 
to  find  e  for  a  specified  value  of  the  mass  M  and  the  mass  eccentricity  etc. 

To  keep  the  development  as  general  as  possible,  nondlmensional  parameters 
will  be  used.  These  are: 

•  The  nondlmensional  unbalance  parameter  a 

3  3 

•  A  nondlmensional  speed  or  mass  parameter  a  »  Muc  /(yRL  ) 

The  significance  of  these  parameters  will  become  clearer  as  analysis  is 
developed.  Ideally  we  will  develop  charts  from  which  e  can  be  directly 
read  when  values  of  a  and  a  are  specified.  Unfortunately,  because  of  the 
nature  of  the  equations  involved,  the  equilibrium  analysis  cannot  proceed 
directly.  Specifically,  to  avoid  any  need  for  iteration,  it  is  necessary 
to  seek,  from  the  equations,  values  of  a  which  correspond  to  given  values  of 
a  and  e.  Of  course,  once  a  table  of  matching  values  of  a,  e  and  a  is 
established  it  matters  little  how  the  values  were  arrived  at,  and  the 
desired  charts  can  be  directly  created  from  the  tables. 

2 

In  this  analysis,  equilibrium  requires  that  the  centrifugal  force  Mu  6c 
acting  on  the  rotor  mass  is  equal  to  the  fluid  film  force  F  acting  on  the 
damper : 

Mu2 8c  -  F  (2-J09) 


So  recalling  the  definition  of  o, 


a 


Fc 


uuRL' 


3  ' 


1 

6  ’ 


(2-11C) 
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and  if,  given  a  and  e,  we  find  6  and  F,  Che  problem  is  solved.  The  total 
fluid-film  force  F  is  related  to  its  components  as  follows: 


F  -Vfr2  +  F,2 


(2-111) 


and,  the  fluid-film  analysis  based  on  short  bearing  assumptions  and  the  it 
or  half -film  treatment  of  film  rupture  gives  the  following  expressions  for 

F_  and  F_  in  terms  of  e  and  the  damper  geometry  parameters  and  frequency 

R  1 

(for  circular  orbit  whirl)  : 


F  C2 

FR  C  .  2...  2.2 

- -  2e  /(1-e  ) 

uRLw 


ft  02 

OT  3 
URL  a) 


ir  e/U(l-e2)3/2) 


(2-112) 


(2-113) 


An  Inspection  of  the  geometry  of  equilibrium  shown  in  Figure  2-6  gives 
relationships  between  fluid  film  components  and  the  load  angle  y: 


—  -  tan  Y 
FR 


(2-114) 


and  the  triangle  whose  sides  are  ac,  Sc  and  ec  gives  the  quadratic  relation¬ 
ship  for  6 : 


222  1/2 
S  ■  e  cos  y  +  (a  -  e  sin  y) 


(2-115) 


With  the  above  equations,  the  mass  parameter  a  can  be  calculated  from  given 
values  of  eccentricity  ratio  e  or  unbalance  parameter  a  using  the  following 
scheme : 


Step 

Number 


Value  Assumed  or  Resulting 
from  Previous  Step _ 


(ft>  fr) 


(ft»  fr) 


(e,  a,  y) 

Fc2  . 

(8,  — rl 

URL  w 


Calculate 


(Fr  V  -V 

tiRL  u 


3 

URL  u 


Using  Equation 


(2-112)  and  (2-113) 


(2-114) 


(2-111) 


(2-115) 


(2-110) 


As  Intended,  this  procedure  yields  sets  of  values  of  a,  a  and  e  from  which 
charts  can  be  prepared.  It  should  be  noted  that,  for  certain  combinations 
of  the  mass  or  speed  parameter  a,  and  the  unbalance  parameter  a,  more  than 
one  equilibrium  value  of  orbit  radius  e  exists  and,  in  fact,  as  many  as 
three  can  exist,  as  explained  at  the  beginning  of  this  section. 


2.6.2  Coupling  to  the  Finite  Length  Damper  Program,  Including  Shaft 
Flexibility  and  Parallel  Springs 

The  squeeze  film  damper  analysis  program  calculates  and  stores  tables  of 

direct  (B  )  and  cross  coupled  (B  )  circular  orbit  damping  coefficients  as 
xx  xy 

functions  of  eccentricity  ratio.  The  coefficients  are  in  nondimens ional 
fora,  specif ically: 


(2-116a) 


yR(L/c)' 


xy  yR(L/c)3 


(2-116b) 


The  interfacing  to  the  program  is  done  by  interpolating  in  the  tables  for 

the  B  and  B  values  corresponding  to  the  eccentricity  value  of  Interest, 
xx  xy 

Nondimens ional  forces  are  then  calculated  using  the  equations: 


£t(s)  ■  B«.s 


(2*-117 a) 


f  (e)  -  B  e 
r  xy 


(2-11 7b) 


The  forces  can  now  be  used  in  an  equilibrium  analysis  which  incorporates 

shaft  (or  bearing)  flexibility  in  series  with  the  squeeze  film.  The  flexible 

element  is  assumed  to  be  light  and  can  be  envisioned  as  a  bearing  which 

carries  a  rigid  rotor  or  a  simple,  Jef f cott-type  rotor  on  stiff  bearings. 

Figure  2-7  is  an  equilibrium  diagram  showing  the  modification  to  the  mass 

center  location  caused  by  a  flexible  element  of  stiffness  R  in  series  with 

s 

the  damper.  If  we  define  f  from  Equation  (2-117)  to  be 


■|tr2  * 


(2-118) 
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Figure  2-7  Damper  Rotor  Equilibrium  for  Case  of  a  Flexible 
Element  in  Series  with  the  Damper 


then  forces  F^,  F^,  and  F  are: 


ft(e)  uiyRL^/C2 

(2-119a) 

frC e)  uyRL^/C2 

(2-119b) 

fCe)  i»)uRL2/C2 

(2-119c) 

From  geometry  (cosine  rule)  we  get, 


/ _ i \ 2  2  r  F  s 2,  cF 

(E  5  £  +  CKC)  +  2  K“C  COSY 

s  s 


(2-120) 


and  from  the  sine  rule  we  have: 


sin  y  * 


(2-121) 


From  equation  (2-119c)  we  get 


F  fwyRL  1 
C2 


(2-] 22) 


and  a  parameter  f  for  shaft  flexibility  can  be  defined  as 


therefore 


.  WR  (U3 

s  K  V 
s 


— — —  •  f  y 
CK  rs 
s 


(2-123) 


(2-124) 


The  effect  of  springs,  such  as  flexures,  in  parallel  with  the  damper  film 
can  also  be  modelled  by  defining  a  parallel  spring  flexibility  parameter 
similar  to  Equation  (2-123) : 


f 


(2-124) 


where  is  the  parallel  spring  stiffness.  The  radial  force  f  must  then 
be  modified  to  account  for  this  additional  force 


f 


r 


fr +f: 


(2-117c) 
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The  equilibrium  computation  then  proceeds  as  follows:  a  value  of  e  is 
assumed  and  values  of  y»  fr»  f£  and  f  are  calculated  using 


,  -1  ,rt, 

Y  -  tan  (y-) 

r 

(2-125) 

and  Equations  (2-117)  and  (2-118).  The  modified  eccentricity  ratio  e'  is  found 
from  Equation  (2-120),  y*  is  calculated  from  Equation  (2-121),  rewritten  as 

,  .  -1  /tsinYv 

y'  -  sin  (  e,  ) 

(2-121a) 

and  force  components  F_  and  F_  are  extracted  from  Equation  (2-119) . 

R  T 

inertia  force  applied  to  the  mass  must  equal  the  fluid  film  force  (as 
modified  by  series  and/or  parallel  stiffnesses)  : 

Now  the 

FT'  -  Ft  cos  (y-y')  -  FR  sin(Y-Y') 

(2-126a) 

•  Fr’  -  Fr  cos  (Y-Y*)  +  ft  sin(Y-Y’) 

(2-126b) 

and  similarly 

f  '  *  f*  cos(Y-Y')  -  f  sin(  Y- Yf) 
t  t  r 

(2-127a) 

fr’  •  fr  COS(Y-Y')  +  ft  sin(Y-Y’) 

(2-127b) 

f*  «  [(fr)2  +  (f,.)2]172 

(2-127c) 

From  Equation  (2-115),  2  values  of  S  (6^  and  0^  can  be  found  to  be: 

2  -  e'  cosY’  +  [(o)2  -  (ef  sinY’)2]172  (2-128) 
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so  that  the  condition  is 


8'  -  r 

(2-131) 

Ue  therefore  get  two  values  for  the  nondlmensional  speed 

or  mass  parameter 

cr: 

f* 

"1“  *1 

(2-132a) 

f ' 

°2  “  $2 

(2-132b) 

Therefore  for  given  values  of  the  displacement  parameter  (or  eccentricity 
ratio)  e  and  mass  eccentricity  parameter  a  we  can  find  equilibrium  solutions, 
if  they  exist,  for  the  speed/mass  parameter  a.  These  values  can  be  plotted 
to  generate  a  "jump  map"  of  the  type  shown  in  Figure  (2-8)  ,  which  is  an 

example  using  short  bearing  analysis,  where  it  was  assumed  that  the  shaft 
was  infinitely  flexible  (JF  ■  Q)  and  the  flexures  were  infinitely  rigid 

CFf  »  -I . 


The  above  equations  have  been  programmed  and  samples  of  the  input  and  output 
are  given  in  Section  III  of  this  report.  An  additional  quantity  calculated 
by  the  program  is  transmlssibility 


(2-133) 


which  is  often  of  interest,  especially  when  force  transmitted  to  pedestals 
is  of  concern. 
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Figure  2-8  Map  of  Displacement  Parameter  Versus  Speed/Mass  Parameter 
for  Different  Values  of  Mass  Eccentricity  Parameter  (a) 


SECTION  III 


DAMPER  ANALYSIS  COMPUTER  PROGRAMS 

The  squeeze  film  damper  analysis  described  in  Section  II  of  this  report  has 
been  implemented  in  the  form  of  a  series  of  computer  programs.  These  pro¬ 
grams  are  linked  by  data  files  to  form  an  interlocking  set  of  tools  for 
generation  and  display  of  damper  performance  data,  and  to  provide  an  inter¬ 
face  with  rotordynamics  system  analysis  programs. 

In  this  section,  the  individual  programs  will  be  described  as  they  would  be 
used  in  the  context  of  a  damper  component  design  analysis  system.  In  addi¬ 
tion,  examples  of  program  input  and  output  will  be  presented  which  exercise 
major  program  options. 

Program  Descriptions 

A  total  of  eight  computer  programs  have  been  written  to  generate  and  process 
squeeze  film  damper  data.  These  programs  are  listed  in  Table  3-1. 


Table  3-1 

Squeeze  Film  Damper  Analysis  Programs 


1.  SQZDMP 

2.  JUMP 

3 .  DIMOUT 

4.  INTERP 

5.  CIRCULAR  ORBIT  GRAPHICS 

6.  POST 

7.  NON-CIRCULAR  ORBIT  GRAPHICS 

8.  INTERFACE 


Figure  3-1  illustrates  the  interaction  and  flow  of  data  between  programs. 
Program  SQZDMP  joins  the  analytical  base  and  must  be  run  Initially  in  all 
cases  to  generate  the  appropriate  data  file(s)  through  the  solution  of  the 
Reynolds  equation.  The  files  generated  can  then  be  accessed  as  shown  in 
Figure  3-1  by  any  of  the  seven  post  processing  routines.  Each  individual 
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Figure  3-1  Squeeze  Film  Damper  Analysis  Program 


program  will  now  be  described  in  functional  grouping  according  to  the  SQZDMP 
option  chosen,  either  dimensionless,  dimensional,  or  mobility. 

a.  Dimensionless  Circular  Orbit  Calculations 

This  option  involves  the  exercise  of  as  many  as  five  programs  and  is 
intended  for  use  as  a  design  tool  wherein  the  effects  on  performance  of  a 
variety  of  damper  parameters  can  be  conveniently  determined.  The  damper 
is  assumed  to  be  executing  a  centered,  circular  orbit  at  a  user-specified 
range  of  eccentricities. 

It  is  important  to  note  that  by  performing  the  initial  calculations  in 
dimensionless  form,  basic  damper  trends  may  be  examined.  For  example, 
damper  clearance,  oil  viscosity  and/or  supply  pressure  may  be  varied 
conveniently  to  determine  their  effects  on  damper  performance  for  a  given 
damper  geometry  without  having  to  repeat  the  damper  fluid-film  solution 
sequence  for  every  such  variation.  Dimensionless  files  also  provide  an 
efficient  way  to  supply  data  for  programs  performing  dynamics  calculations 
(such  as  nonlinear  steady-state  damped  response) . 

SQZDMP,  regardless  of  the  option  chosen,  requires  as  input  a  definition  of 
damper  geometry,  which  consists  of  the  feeding  and  sealing  arrangement  and 
the  finite  difference  mesh  used  for  calculations,  the  leakage  and  flow 
coefficients  and  the  necessary  program  control  parameters.  The  dimension¬ 
less  circular  orbit  option  in  particular  requires  the  L/D  ratio,  nondimen- 
sional  supply  and  cavitation  pressures  and  a  series  (up  to  12)  of  values 
for  journal  center  eccentricity  and  angular  location  with  respect  to  the 
coordinate  axes  (in  order,  for  example,  to  locate  the  journal  center  with 
respect  to  a  feed  or  drain  hole) .  The  result  of  the  calculation  is  a  table 
of  dimensionless  circular  orbit  coefficients  which  is  also  suitable  for  use 
in  a  design  chart.  This  information  includes: 

e  dimensionless  damping 
e  dimensionless  inertia 
e  dimensionless  flow  ln/out 
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•  dimensionless  pressure  distribution  (optional) 

and  is  put  out  in  the  form  of  printed  tables  or  a  data  file  to  be  accessed 
automatically  by  other  programs. 

CIRCULAR  ORBIT  GRAPHICS  generates  plots  of  the  nondimens ional  SQZDMP  - 
generated  damping  and  inertia  coefficients  as  functions  of  eccentricity 
ratio.  This  allows  for  convenient  checking  of  the  calculations  and  for 
their  use  as  design  charts.  Input  quantities  are  the  tables  of  dimensionless 
coefficients  produced  by  SQZDMP  and  the  appropriate  plotting  control  variables. 
Program  output  consists  of  plots  of  damping  or  inertia  coefficients  (direct 
or  cross-coupled)  as  functions  of  eccentricity  ratio. 

DIMOUT  uses  as  input  the  nondimens ional  circular  orbit  coefficient  tables 
generated  by  SQZDMP  and  dimensionalizes  them  using  a  set  of  user-input 
damper  dimensions.  Thus  a  single  nondimens ional  file  may  be  used  to 
represent  a  range  of  dampers  (of  similar  configuration  but  different  sizes, 
etc.)  during  a  design  exercise  without  having  to  rerun  the  analysis  program 
repeatedly. 

INTERP  accepts  the  dimensionless  data  generated  by  SQZDMP  along  with  supply 
pressures  or  eccentricity  ratios  and  provides  as  output  dimensionless 
circular  orbit  coefficient  tables.  The  program  uses  an  interpolation 
scheme  so  that  coefficient  values  may  be  provided  which  are  in  between 
those  in  the  input  tables.  INTERP  can  also  be  used  to  couple  a  damper  model 
with  a  rotordynamics  computer  program.  For  example,  a  nonlinear  steady-state 
response  program  can  extract  successive  values  of  damping  coefficients  as  it 
iterates  on  response  amplitude  by  using  INTERP  to  access  the  appropriate 
nondlmenslonal  data  files  and  interpolate  to  find  the  correct  coefficients 
for  a  given  damper  circular  orbit  eccentricity. 

JUMP  uses  as  input  a  dimensionless  circular  orbit  coefficient  file  produced 
by  SQZDMP  along  with  parameters  which  describe  damper  displacement,  mass 
eccentricity  and  shaft  flexibility  and  parallel  spring  stiffness  if  present. 

The  program  output  is  a  series  of  tables  which  relate  the  damper  orbit 
circle  radius  and  mass  eccentricity  with  a  third  parameter  representing 
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nondimens ional  speed  or  mass.  These  results  are  plotted  manually  to  form 
"jump  maps"  which  reveal  regions  where  the  damper  may  be  bistable,  that  is 
have  two  possible  orbit  radii  at  which  stable  operation  is  possible. 

b.  Dimensional  Calculations 

Calculations  for  a  squeeze  film  damper  can  also  be  performed  based  on  a 
specific  set  of  user  supplied  damper  dimensions  and  operating  conditions. 
Data  file  output  is  not  provided  for  this  option  since  it  is  designed  for 
evaluations  of  individual  dampers  rather  than  to  study  variations  of  a 
design.  With  this  calculation  the  damper  journal  position  and  velocity  are 
specified  for  each  eccentricity  by  the  user.  Therefore  circular  orbit  con¬ 
ditions  can  be  simulated  as  well  as  individual  points  on  an  arbitrary 
noncircular  orbit ,  for  which  instantaneous  values  of  damping  and  Inertia 
coefficients  (valid  for  those  points  only)  are  provided. 


SQZDMP  (Dimensional  Option)  is  the  only  program  which  is  run  under  this 
option.  In  addition  to  the  standard  set  of  input  data  described  previously, 
it  requires  the  damper  dimensions,  a  reference  operating  frequency,  viscosity 
and  density  of  the  damper  fluid,  dimensional  supply  and  cavitation  pressures 
and  the  location  and  velocity  of  the  damper  journal.  The  output  includes 
Instantaneous  values  of  force  and  local  damping  coefficient,  circular  orbit 
damping  and  inertia  coefficients,  power  loss  and  pressure  distributions  and 
is  available,  as  noted  above,  as  a  printout  only. 

c.  Mobility  Calculations 

Up  to  this  point,  calculations  have  been  described  which  are  either  steady- 
state  (circular  orbit)  or  quasi-static  (calculations  for  a  given  journal 
position  and  velocity  vector).  However,  it  is  also  necessary  to  be  able  to 
calculate  damper  behavior  in  response  to  an  arbitrary  load  history.  This 
must  be  a  time-transient  calculation  since  the  forces  on  the  damper  can 
Include  both  steady  and  time-varying  components  and  the  inherent  nonlinearity 
of  the  damper  causes  it's  response  to  depend  heavily  on  the  damper  journal’s 


position  and  velocity  vector.  To  increase  the  efficiency  of  the  calculation 
the  nobility  method,  described  in  Section  2.5,  is  used.  This  method  avoids 
the  necessity  of  re-solving  the  Reynolds  equations  for  every  time  step  by 
basing  damper  response  on  nondimens ional  pressure  tables  which  have  been 
precalculated  by  SQZDMP.  Thus,  a  laborious  equation  solution  scheme  has 
been  replaced  by  a  much  simpler  table  -  lookup  procedure. 

A  total  of  four  programs  may  be  exercised  under  this  option.  Post-processing 
involves  both  damper  response  analysis  and  also  linkup  of  the  pressure  tables 
to  an  external  rotordynamics  program  (user-supplied)  for  system  time- 
transient  calculations. 

SQZDMP  (Mobility  Option),  which  performs  dimensional  calculations  only, 
requires  the  same  input  quantities  as  the  dimensional  program  option,  except 
that  the  journal  center  coordinates  are  read  in  as  a  set  (12  maximum)  to 
produce  the  pressure  tables  corresponding  the  journal  eccentricity  points 
for  unit  X  and  Y  velocities.  The  pressure  tables  may  be  used  in  two  ways. 
First,  they  can  be  used  in  the  calculation  of  damper  force  during  time 
transient  analyses  using  a  rotordynamics  program,  and  secondly  they  can  be 
used  for  rigid-rotor  time  transient  analysis  using  a  built-in  mobility- 
based  calculation  scheme.  For  the  latter  calculation,  a  load  history  must 
be  specified  along  with  a  convergence  error  and  maximum  iteration  number. 

The  program  will  then  use  the  mobility  method  to  compute,  in  printed  table 
and  data  file  form,  a  response  time  history  consisting  of  force,  displace¬ 
ment  and  velocity  as  functions  of  time,  based  on  the  pressure  tables 
generated,  and  input  parameters  describing  the  applied  load. 

POST  calculates  the  following: 

1.  Effective  stiffness  and  damping  based  on  the 
time-transient,  noncircular  orbit  calculation 
results. 

2.  Frequency  spectra  of  relevant  quantities  (e.g. 
displacement)  to  determine  synchronous  and 
non-synchronous  content . 


Results  are  presented  in  tabular  form. 


NON-CIRCULAR  ORBIT  GRAPHICS  displays  plots  of  the  time-transient  mobility 
method  calculation.  It  can  display: 


•  hysteresis  loops 

•  displacement  or  force  as  a  function  of  time 

•  journal  orbits  (X  vs  Y  amplitude) 

•  frequency  spectra  of  calculated  waveforms. 

INTERFACE  couples  the  pressure  tables  generated  by  the  mobility  option  of 
SQZDMP  with  an  external  rotordynamics  program  (e.g.  a  time-transient  analysis 
program).  Journal  velocity  and  displacement  information  is  given  to  INTERFACE, 
which  then  returns  the  corresponding  damping  and  inertia  coefficients. 

2.  Sample  Input/Output 

Figure  3-2  is  an  expanded  version  of  Figure  3-1  and  shows  in  more  detail 
the  interaction  of  the  damper  analysis  programs  and  their  individual  input 
requirements  and  output  quantities.  The  programs  themselves  are  shown  in 
heavily  outlined  boxes,  required  input  quantities  in  lightly  outlined  boxes 
and  program  output  quantities  are  shown  in  parallelograms.  Solid  lines 
denote  data  flow  through  a  given  program  and  broken  lines  show  data  trans¬ 
fer  by  file. 

To  illustrate  program  Input/output,  an  example  case  for  an  end  groove 
feed/end  groove  drain  damper  was  prepared.  This  damper  geometry  is  shown 
schematically  in  Figure  3-3  along  with  the  relevant  dimensions  and  operating 
parameters.  The  balance  of  Section  3.2  shows  calculations  made  for  this 
damper  using  dimensionless,  dimensional  and  mobility  options. 

a.  Dimensional  Calculations 

The  damper  configuration  described  in  Figure  3-3  was  first  analyzed  using 
the  dimensional  input  option  to  SQZDMP.  This  is  highlighted  in  Figure  3-4, 
which  has  been  extracted  from  Figure  3-2  and  illustrates  the  flow  of 
information  for  this  option. 


END  GROOVE  FEED /DRAIN  DAMPER 


"Global"  Seal 


PARAMETERS 

Length  * 

Diameter  * 

Radial  Clearance  ■ 

Fluid  Viscosity  ■ 

Fluid  Density  * 

Frequency  ■ 

Feed  Groove  Pressure  ■ 

Cavitation  Pressure  ■ 

Journal  Coordinates  ■ 

Journal  Velocities  ■ 

Seal  Flow  Coefficient  * 


1.0  in 
5.0  in 
.008  in 

1.0  x  10  ®  lbs-sec/in^ 
7.9  x  10  lbs-sec^/in^ 
133.33  Hz  (8,000  RPM) 

20  pslg 
14.7  psig 

X  -  0.004  in 

Y  -  0 

V  -  0 

V*  “  3.35  in/sec 

2.0 


Figure  3-3  Sample  Problem  for  Exercising  Damper  Software 
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SQZDMP 


Lrcular  Orbit  Option 


Dimensional  Option 
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Output  from  the  dimensional  calculation  includes: 

•  Damper ■ forces 

e  Instantaneous  damping  coefficients  (for  the 
specified  journal  position  and  velocity) 

e  Instantaneous  power  loss 

•  Tondl  coefficients 

•  Circular  orbit  damping  coefficients 

•  Circular  orbit  inertia  coefficients 

•  Flows  into /out  of  damper 

•  Pressure  distribution. 

Table  3-2  shows  the  complete  input  and  output  for  a  dimensional  calculation 
using  the  damper  geometry  and  operating  parameters  from  Figure  3-3  as  input 
quantities. 

b.  Dimensionless  Calculations 

The  next  analysis  to  be  performed  exercises  the  dimensionless  option  of 
SQZDMP  and  key  post-processing  programs.  The  flow  of  information  for  this 
option  is  shown  in  Figure  3-5  (again,  a  subset  of  Figure  3-2).  The  damper 
to  be  analyzed  is,  as  before,  that  described  in  Figure  3-3  but  because  of 
the  nature  of  the  calculations  nondimenslonalization  of  some  input  parameters 
was  required  and  a  series  of  eccentricity  ratios  was  specified  (as  opposed 
to  a  single  eccentricity  the  dimensional  case).  The  new  parameters  are: 

L/D  ratio  ■  0.2 

Supply  pressure  ■  0.02122 

Cavitation  pressure  ■  -0.0156 

Journal  center  coordinates 

(12  values)  X  *  .00001, .1, . 2, . 3, .4, . 5 , .6 , . 7 , .8 , .9 , .95, .99 

Y  •  0.0  (all  values) 

From  this  input  data,  SQZDMP  generates  an  "echo  print"  of  the  input  data,  a 
series  of  12  tables  (one  for  each  eccentricity)  of  nondimensional  pressure 
distributions  and  a  summary  table  of  dimensionless  circular  orbit  coeffi¬ 
cients.  Table  3-3a  shows  the  input  data  printout  and  a  representative 


Program  Output  for  Dimensional  SQZDMP  Run 
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Input  Data  Listing  for  Dimensionless  SQZDMP  Run 


pressure  distribution  (for  an  eccentricity  ratio  of  0.5)  is  shown  in 

Table  3-3b,  while  Table  3-3c  gives  the  summary  of  circular  orbit  coefficients 

for  this  case. 


Plots  of  dimensionless  damping  and  inertia  are  available  through  the  use 
of  the  CIRCULAR  ORBIT  GRAPHICS  program.  Figures  3-6a  and  3-6b  are  plots  of 
direct  and  cross-coupled  damping  coefficients  for  the  output  printed  in 
Table  3-3c. 

The  dimensionaless  data  generated  by  SQZDMP  can  also  be  dimensionalized 
using  program  DIMOUT.  Table  3-4  shows  a  typical  DIMOUT  summary  table  for 
the  example  damper. 

The  results  of  running  the  JUMP  program,  for  a  bi-stable  orbit  or  "jump" 
analysis,  are  given  in  Table  3-5.  For  this  case  both  shaft  and  flexure 
flexibility  were  assumed  to  be  zero  and  the  resulting  output  gives  values 
for  the  speed/mass  parameter  "SIGMA"  in  terms  of  the  displacement  parameter 
(or  eccentricity  ratio)  "EPS"  and  the  mass  eccentricity  parameter  "ALPHA" 
along  with  intermediate  results  and  the  transmissibility  "T". 

c.  Mobility  Calculations 

The  final  analysis  to  be  performed  uses  the  Mobility  option  to  SQZDMP  and 
its  key  post  processing  programs  as  described  in  Figure  3-7.  Since  the 
Mobility  calculation  is  dimensional,  the  damper  test  case  described  in 
Figure  3-3  was  run  for  a  new  set  of  journal  center  coordinates  (to  generate 
the  required  pressure  tables)  of: 

Journal  center  X  ■  .00001,  .001,  .004,  .006,  .0079  in 

eccentricity  Y  -  0.0,  0.0,  0.0,  0.0,  0.0 

For  the  mobility  calculation,  the  Journal  was  set  initially  at  X  ■  .0001  in., 
y  «  -.002  in.  (i.e.  slightly  off-center).  A  steady  state  force  of  25  lbs. 
was  applied  in  the  X  direction  and  a  rotating  unbalance  of  8.5  x  10  ^  lbs- 
sec^  at  260°  was  also  applied.-  Table  3-6a  shows  the  input  data  for  this 
calculation  and  Table  3-6b  is  the  first  page  of  the  calculation  output 
showing  the  stepwise  time  integration. 
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Pressure  Distribution  for  Dimensionless  SQZDMP  Run 
(Eccentricity  Ratio  =  0.5) 


TWcftiWr  ft  A1  ft  ft  ftfrft  ft  ftftftft  jHfiHIp  ft  ft^rftft  ft  ft  ft  ft  ft  -A  A  ft  ft  'A'  ftftftft 

*  SQUEEZE  FILM  DAMPER  PERFORMANCE  TABLE  *  ********** 

*  CIRCULAR  ORBIT  DATA  *  *01/07/8 2* 

*  MTI  PROGRAM  SQZDAMP  (PN812)  *  ********** 

********************  *  ** *  ******  ***** ****** 

CASE  2  -  SEALED  DAMPER  -  DIMENSIONLESS  OPTION  1/7/82  F.GILLHAM 


ftft  ft  ftftftft  NON-DIMENSIONAL  DATA  ftftftft  ftft  ftft 

L/D  RATIO  .  0.2000 

SUPPLY  PRESSURE  PARAMETER  (PSP/ (12W*MU*(R/C)**2) ) , .  2.1220D-02 

CAVITATION  PARAMETER  (PCV/ (12W*MU*(R/C)**2) )  .  -1.5600D-02 

LOCAL  SEAL  PARAMETER  1  (CL1*12MU*R/C**3) .  0.0 

LOCAL  SEAL  PARAMETER  2  (CL2*12MU*R/C**3) .  0.0 

OIL  SUPPLY .  CIRCUMFERENTIAL  GROOVE 

DAMPER  IS  .  .  ASYMMETRIC 


ftftftftft  ftft  ftft  CIRCULAR  ORBIT  COEFFICIENTS  ********* 


ORBIT  |  *****  DAMPING  *****  |  *****  INERTIA  *****  |  ****  FLOW  *** 


RADIUS/C  |  B/(MU*L*(R/C)**3)  1  INERT. / (RH0*L*R**3/C) |  Q/(W*C*R**2) 


|  DIRECT 

I (BXX=BYY) 

I 

X-COUPLED 

(BXY=-BYX) 

|  DIRECT 
|  (CXX=CYY) 

1 

X-COUPLED 

(CXY=-CYX) 

|  (SEE  NOTE  1)  | 

|  IN  OUT  | 

1 _  1 

0.000010 

4 . 886D-01 

-4.334D-07 

4.886D-02 

4.153D-14 

3.8D-02 

3.8D-02 

0.100000 

4.957D-01 

-4.306D-11 

4.92 3D-02 

3 . 780D-14 

3.8D-02 

3.8D-02 

0.200000 

5 . 182D-01 

-2 . 149D-11 

5.038D-02 

3.470D-14 

3.8D-02 

3.8D-02 

0.300000 

5 . 597D-01 

-1.435D-11 

5 . 244D-02 

3.206D-14 

3.8D-02 

3.8D-02 

0.400000 

6 . 280D-01 

-1.080D-11 

5.568D-02 

2.980D-14 

3.8D-02 

3.8D-02 

0.500000 

7 . 386D-01 

-8.674D-12 

6.06  ID-02 

2. 784D-14 

3.9D-02 

3.9D-02 

0.600000 

9 . 084D-01 

2.680D-02 

5 . 884D-02 

5 . 475D-03 

3.9D-02 

3.9D-02 

0.700000 

1. 103D+00 

3.174D-01 

5 . 633D-02 

1 . 337D-02 

5.0D-02 

2.9D-02 

0.800000 

1 . 482D+00 

1 . 184D+00 

7.015D-02 

1 . 908D-02 

2.9D-02 

5.2D-02 

0.900000 

2 . 863D+00 

4 . 852D+00 

1. 111D-01 

3 . 079D-02 

1.3D-02 

6.9D-02 

0.950000 

5 . 844D+00 

1.413D+01 

1. 862D-01 

4.936D-02 

-5.0D-03 

8.7D-02 

0.990000 

1 . 737D+01 

5 . 389D+01 

8.295D-01 

1 . 632D-01 

-1.8D-02 

9.9D-02 

NOTE  1:  DUE  TO  CAVITATION  CONDITION  ASSUMPTLONS, 

PREDICTED  FLOW  IN  DOES  NOT  NECESSARILY  EQUAL  FLOW  OUT 

ABBREVIATION  KEY  :  W  *  ROTATIONAL  SPEED  R  =  DAMPER  RADIUS 

MU  *  FLUID  VISCOSITY  C  *  DAMPER  CLEARANCE 

RHO  *  MASS  DENSITY  L  *  DAMPER  LENGTH 


Table  3-3c 


Summary  Output  for  Dimensionless  SQZDMP  Run 


CIRCULAR  ORBIT  COEFFICIENT  -  BXX  -  BVV 


SEALED  DAMPER  -  CIRCULAR  ORBIT  COEFFICIENTS  -  BXV 


*  DIMENSIONAL  OUTPUT  FOR  SQUEEZE  FILM  DAMPER  PERFORMANCE 

*  CIRCULAR  ORBIT  OATA 


CASE  2  -  SEALED  DAMPER  -  DIMENSIONLESS  OPTION  1/7/82  F.GILLHAM 
01/07/82 


DIAMETER 

LENGTH 

CLEARANCE 

VISCOSITY 

FREQUENCY 

DENSITY 

( 1 NCH ) 

(INCH) 

( INCH) 

(LBF-SEC 
/ 1 N**2 ) 

(HERTZ) 

( L8F-SEC**2 
/ 1 N**4 ) 

5 . 0000 E *00 

1.0000E+00 

8.0000E-03 

1 . 0000E-06 

1 .3333E+02 

7. 9000E-05 

L/0  RATIO  .  0.2000 

SUPPLY  PRESSURE  <LBF/IN*#2)  2.0833E+01 

CAVITATION  PRESSURE  (LBF/IN**2)  -1.5315E+01 

GLOBAL  SEAL  PARAMETER  1  ( IN**3/SEC/PSI )  .  8.5333E-02 

GLOBAL  SEAL  PARAMETER  2  ( I N**3/SEC/ PS  I )  .  8.5333E-02 

LOCAL  SEAL  PARAMETER  1  ( I N**3/SEC/ ( PSI - I N ) ) .  0.0 

LOCAL  SEAL  PARAMETER  2  ( I N**3/S£C/( PS  I - 1 N ) ) .  0.0 


OIL  SUPPLY . CIRCUMFERENTIAL  GROOVE 

OAMPER  IS  . ASYMMETRIC 


CIRCULAR  ORBIT  COEFFICIENTS  ***** 


ORBIT 

DIRECT 

X-COUPLED 

DIRECT 

X-COUPLED 

FLOW 

FLOW 

RAO 1  US 

DAMP  1 NG 

OAMP 1 NG 

INERTIA 

INERTIA 

IN 

OUT 

INCH 

LBF-SEC 

LBF-SEC 

LBF-SEC**2 

LBF-SEC**2 

1  N**3 

1  N**3 

/IN 

/IN 

/IN 

/IN 

/SEC 

/SEC 

0.0 

1.491E+01 

-1.323E-05 

7.539E-03 

6.408E-15 

1 . 577E+00 

1 . 577E+00 

0 . 800E-03 

1.513E+01 

-1 . 314E-09 

7 . 596E-03 

5.832E-15 

1 . 580E+00 

1 . 580E+00 

0.  160E-02 

1 . 581E+01 

-6.559E-10 

7.773E-03 

5.353E-15 

1 . 587E+00 

1 . 587E+00 

0.240E-02 

1.708E+01 

-4.379E-10 

8.092E-03 

4.947E-15 

1 . 598E+00 

1.598E+00 

0. 320E-02 

1.916E+01 

-3.295E-10 

8.592E-03 

4.598E-15 

1.612E+00 

1.612E+00 

0 . 400E-02 

2.254E+01 

-2. 647E-10 

9.351E-03 

4.295E-15 

1 . 6Z7E+00 

1 . 627E+00 

0 . 480E-02 

2.772E+01 

8. 180E-01 

9.079E-03 

8.448E-04 

1.642E+00 

1 . 642E+00 

0.560E-02 

3.367E+01 

9-68 6E+00 

8. 69 IE-03 

2.064E-03 

2. T12E+00 

1.231 £*00 

0 . 640E-02 

4. 523E+01 

3.613E+01 

1 .082E-02 

2.943E-03 

1 . 197E+00 

2. 173E+00 

0. 720E-02 

8.736E+01 

1.481E+02 

1 .714E-02 

4.750E-03 

5.409E-01 

2.895E+00 

0.760E-02 

1 . 783E+02 

4.313E+02 

2. 873E-02 

7. 617E-03 

-2.099E -01 

3.647E+00 

0 . 792E-02 

5. 302E+02 

1 . 644E+03 

1 . 28QE-Q1 

2.518E-02 

-7. 337E-01 

4. 167E+00 

PSI(OEG)  *  H0LE1  TO  ECCENTRICITY  VECTOR  IN  WHIRL  DIRECTION 


Table  3-4 

Dimensionalized  Circular  Orblc  Coefficients 
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—MT  fnak/aaol/spoontry  itt latcloo 
HMh  Definition 
Lanka*#  aak  floe  Coef fit lent* 

Control  Htn«m 


i 

Miktlity  Option  (dlMUlantl) 


Stank* rk  Option 


look  hlotorp 
bonttmea  error 
Nu  no.  <9!  Itera¬ 
tions 


Pressure  Table 
Option  BO* 


INPUT  DATA  LISTING 
PROGRAM  PN812 


02/16/82  15:21:25 


. . 

CAR 01  TITLE 

CASE  1  *  SEALED  DAMPER  :  MOBILITY  CALCULATION  AT  8000.  RPM 


CAR02  INTEGER  CONTROLS 
M  N  NOIM  INERT  NXY  INP  (UNITS  ICIRC  MCONST 
24  9  0  0  -5  1  0  0.1 


CAR03B  BEARING  DIMENSIONS 

01 AM  BLG  CLRC  VI  SC  OENS  FREQ 
5.0  1.0  .008  1 . 0000E-6  7.90E-5  133.33 


CARD4  BEARING  DESCRIPTION  OPTIONS 

MSL  MFD  MHF  NFS 
10  10 


CAR09  (SUBMIT  IFMSL-1) 

CGI  CG2 

2.0  2.0 


CARO  10  LUBRICANT  INFORMATION 

PSP  PCV 
20.  -14.70 


CARO  1 1 A 

(8C10.3) 

X-COORDINATE 

OF  JOURNAL  CENTER 

TXY( 1 ) 
.00001  . 

.001  .004  .006  .0079 

CARO  1 1B 

(8E10.3) 

Y-COORDINATE 

OF  JOURNAL  CENTER 

TYZ(I) 

.0 

o 

b 

o 

o 

CARO  13 

(1615) 

USE  FOR: 

MOBILITY  ONLY 

NTST 

MX IT  ITOP  NSTEP 

1  DIG 

240 

25  3  0 

0 

CARO  14 

(8E10.3) 

USE  FOR: 

MOBILITY  ONLY 

XZ  YZ  XV  YV  CVER 

.0001  -.002  0.  0.  .0001 


CARD  17  TIME-HISTORY  DATA 

TIME  DELT  WXINIT  WYINIT  UNBAL  PHI 
0.  •  .0001  25.  0.  8.5E-5  270. 

CASE  1  -  SEALED  DAMPER  :  MOBILITY  CALCULATION  AT  8000.  RPM 
02/16/82 

24  9  0  0  -5  1  0  1 

0.2000+00  0.5000+01  0.1000+01  0. 8000-02  0. 1000-05  0.1330+03  0.7900-04 
10  10 
0.0  0.0 


Table  3~6a  Input  Data  for  Mobility  Calculation 
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Table  3-6b  Output  (First  Page)  fro«  Mobility  Calculation 


Plots  of  this  output  obtained  by  running  the  NON-CIRCULAR  ORBIT  GRAPHICS 
program  are  shown  in  Figures  3-8.  Figure  3-8a  shows  the  damper  journal 
orbit  within  its  housing  for  several  cycles  of  motion.  Figures  3-8b  and 
3-8c  are  plots  of  displacement  and  force  along  the  Y-axis  as  functions  of 
time,  while  Figure  3-8d  is  a  hysteresis  plot  along  the  Y-axis  (force  as  a 
function  of  displacement) . 


CASE  1  -  SEALED  DAMPER  t  MOBILITY  AT  8000.  RPM 


SEALED  DAMPER  «  MOBILITY  AT  8000.  RPM 


SECTION  IV 


INTEGRATION  OF  COMPONENT  ANALYSIS  TOOLS 

An  essential  feature  of  the  application  of  a  damper  component  analysis  tool 
in  the  design  of  engine  dampers  is  its  integration  with  rotordynamics 
analysis  software.  Therefore,  a  key  aspect  of  this  program  has  been  to 
develop  a  coupled  rotor  system/damper  component  analysis  capability.  This 
capability  was  developed  through  the  combination  of  an  existing  rotor- 
dynamic  analysis  program  by  Detroit  Diesel  Allison  (DDA)  with  the  damper 
component  analysis  software  developed  by  MTI. 

This  and  the  following  section  describe  the  interface  between  the  system 
and  component  programs.  Descriptions  of  the  solution  algorithms  used  in 
the  rotor  system  analysis  program  are  presented  to  the  extent  that  treatment 
of  the  interface  logic  can  be  readily  understood.  A  sample  problem  is 
shown  to  illustrate  the  input  and  output  of  the  analysis  program.  It  is 
not  the  intent  of  this  report  to  discuss  the  merits  of  each  analysis  and 
its  applicability  to  coupled  design.  Rather,  the  intent  is  to  provide  the 
reader  a  description  of  the  interface  with  the  coupled  calculations. 


1.  Overall  Integration  Philosophy 

DDA's  rotor  dynamic  analysis  program,  BB60,  is  coupled  to  the  MTI  damper 
analysis  as  shown  in  Figure  4-1.  The  overall  philosophy  is  as  follows: 

1.  Calculate  the  damper  component  output  in  non-dimensional 
form  and  store  in  files.  Up  to  4  dampers  may  be  calcu¬ 
lated  with  2  files  created  for  each.  The  1st  file  contains 
a  table  of  damping  coefficients  as  a  function  of  eccentri¬ 
city  ratio  e/c.  This  file  is  used  for  single  plane  calcu¬ 
lations  where  circular  orbits  are  expected.  The  2nd  file 
contains  a  table  of  velocity  and  supply  pressure  distri¬ 
butions  which  are  also  a  function  of  e/c.  This  file  is 
used  for  2  or  3  plane  calculations  where  non-circular 
orbits  are  expected. 

2.  Generate  dimensional  damper  forces.  Thu  rotor  dynamic 
analysis,  BB60,  uses  the  damper  component  tables  to 
generate  modal  forces.  Input  to  BB60  includes  the  rotor 
mass/elastic  data  and  dimensional  damper  information. 

Damper  size,  oil  properties,  and  oil  film  radial  clearance 
are  used  to  calculate  a  dimensional  force.  The  damper  file 
number,  damper  location  and  system  connection  information 
is  also  given  in  BB60  input  to  locate  the  appropriate  table 
and  to  apply  the  modal  forces  at  the  proper  location. 
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BB60 

EXECUTIVE 

PROGRAM 


-  NONLINEAR  RESPONSE 

OR 

-  TIME  TRANSIENT 
(CIRCULAR  ORBITS  EXPECTED) 


2.  OR  5  PLANE  ANA1  YSTS 

-  TIME  TRANSIENT 

(NON-CIRCULAR  ORBITS 
EXPECTED) 


DAMPER  COMPONENT 
OUTPUT 


COUPLED 

OUTPUT 


♦ 


Figure  4-1  Overall  Integration  Philosopny 
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3.  Calculate  rotor /damper  response.  The  damper  motion  Is 
unknown  in  general  and  must  be  determined  by  an  Itera¬ 
tion  scheme  so  that  the  damper  calculated  motion  matches 
that  used  to  determine  its  force.  For  the  time  transient 
analysis  the  force  is  determined  from  the  previous  time 
step  damper  motion.  Details  are  discussed  in  the  rotor/ 
damper  interface  section  of  this  report.  In  all  calcu¬ 
lations  the  damper  forces  are  interpolated  from  the 
circular  orbit  coefficients  for  the  single  plane  analysis 
and  from  the  pressure  tables  for  the  2  or  3  plane  analysis. 


The  damper  component  analysis  and  the  BB60  analysis  may  be  executed  separately 
or  together.  The  damper  component  analysis  may  be  executed  in  its  dimensional 
or  non-dimensional  form.  Coupling  the  damper  to  BB60  in  its  non-dimensional 
form  allows  a  one  time  calculation  of  the  damper  model  and  saves  the  dimen¬ 
sional  quantities  for  BB60  input  where  they  may  be  used  parametrically. 

These  inputs  include: 

1.  Damper  location  in  terms  of  engine  station,  system 
connection,  and  damper  file  number. 

2.  Damper  diameter  and  radial  damper  clearance. 

3.  Oil  viscosity  and  density. 

4.  Supply  pressure. 

The  goal  of  this  integration  philosophy  was  to  provide  for  the  most  direct 
and  simplest  coding.  The  input  of  dimensional  damper  data  in  BB60  allows 
for  rapid  parametric  studies  that  would  be  required  in  any  design  analysis. 

A.  description  of  the  creation  of  BB60  damper  execution  for  the  non-linear 
squeeze  film  dampers  is  given  in  Figure  4-2.  It  shows  how  information 
concerning  the  squeeze  film  dampers  enters  the  BB60  analysis. 


2.  Rotor  Dynamic  Analysis  Tool 

A  general  description  of  the  DDA  rotor  dynamic  analysis  program  (BB60)  is 
presented  to  establish  the  computing  capabilities  and  briefly  describe  the 
solution  technique  used. 
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Figure  u-i  Creation  of  BB60  damper  files. 
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FILE  NUMBER 
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INDEX  NUMBER 

BXY 

BXX 

CXY  -I  INERTIA 
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ECCENTRICITY  RATIO 


NUMBER  OF  AXIAL  STATIONS 
NUMBER  OF  CIRCUMFERENTIAL 
STATIONS 
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-  SUPPLY  PRESSURE 

-  ECCENTRICITY  RATIO 


Figure  4-2  ^ continued) .  Creation  of  8860  damper  files. 
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Program  Capabilities 


DDA  haa  developed  a  computer  program  for  the  analysis  of  coupled  rotor/case 
dynamic  systems.  This  program  is  routinely  used  to  compute  the  critical 
speeds  and  dynamic  response  characteristics  of  an  engine  and/or  component 
rig  with  as  many  as  eight  major  systems  (i.e.,  rotors  and  casings).  Beam 
element  representations  for  the  cone,  cylinder,  and  beam  structures  are 
normally  used,  although  provisions  are  made  for  externally  supplying  shell 
and  platelike  flexibilities  to  the  analysis.  Gyroscopic  effects  are  included 
for  multiple  rotor  systems  with  co-rotational  and  counter-rotational  shafts. 
Non-linear  effects  can  be  included  in  the  dynamic  response  calculations. 

The  program  has  been  structured  to  be  very  user  oriented.  Inputs  are 
simply  generated  from  layout  or  detail  drawings  in  the  form  of  inside  and 
outside  diameters  as  a  function  of  axial  position.  Masses  and  stiffnesses 
of  the  structural  elements  are  computed  internally.  Output  data  are  in  the 
form  of  tabular  frequencies,  mode  shapes,  and  responses  as  well  as  machine¬ 
generated  plots.  The  output  can  also  be  routed  to  microfiche  if  desired. 

The  program  is  written  in  Fortran  IV  and  machine  language  for  the  IBM  370 
system. 

The  basic  program  structure  uses  lumped  parameter  analyses  for  rotor  and 
casing  structures.  Masses  and  Inertias  can  be  computed  internally  or  input 
separately.  A  flexibility  matrix  between  connected  masses  is  computed,  using 
beam  and/or  tapered  cone  equations.  Local  flexibilities  are  accounted  for 
separately;  that  is,  computed  or  tested  externally  and  input  as  a  flexi¬ 
bility  matrix.  Nonbeamlike  structures  are  accommodated  by  the  input  of  a 
modal  description  of  the  ID  modes  of  that  structure  which  have  been 
determined  from  a  finite  element  model.  This  means  that  only  the  center-line 
deflection  is  considered.  Stiffness  coupling  elements  between  substructures 
may  be  either  linear  or  nonlinear. 

The  undamped  critical  speeds  of  linear  systems  and  the  damped  linear  and/or 
nonlinear  steady-state  response  to  rotor  unbalances  are  computed  by  the  rotor 
analysis  program.  Transient  solutions  are  determined  for  various  ground  and 
engine  shocks  as  well  as  response  to  a  suddenly  applied  unbalance,  such  as 
blade  loss. 
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Ocher  capabilities  of  che  analysis  cool  include: 


•  Damped  naCural  frequencies 

•  Resonances  at  constant  rotor  speeds 

•  Whirl  flutter  stability 

•  Loads  due  to  misalignment  of  redundant  bearing 
arrangements 

•  System  deflections  resulting  from  discrete  forces, 
moments,  and  loads  due  to  gravity. 

The  rotor  dynamics  program  is  also  used  to  compute  stall-induced  rotor/case 
response  characteristics.  The  case-to-rotor  relative  displacement  is  cal¬ 
culated  for  a  force  on  the  rotor  and  equal  but  opposite  force  on  the  case, 
where  the  force  rotates  at  a  speed  subsynchronous  to  the  rotor  speed.  Again, 
the  normal  case  definition  is  beamlike  but  ID  shell  modes  can  be  accommodated 
with  a  modal  input  description. 

Solution  Technique 

The  prediction  of  rotor  dynamic  characteristics  is  performed  in  BB60  using 
a  modal  synthesis  technique  Qj,  aT]^  as  shown  in  Figure  4-3. 

Step  1  -  The  input  data  describing  the  simulation  are  read.  These  data 
Include  the  description  for  each  substructure  (rotors  and  case)  and  the 
interconnecting  mount  characteristics  (stiffness  and  damping) .  Calculation 
mode  selections  are  read. 

Step  2  -  Mass  and  flexibility  matrices  are  computed  for  each  substructure 
and,  assuming  the  end  mass  stations  to  be  simply  supported  (pinned),  the 
uncoupled  modes  (natural  frequencies  and  associated  mode  shapes)  are  com¬ 
puted  without  gyroscopic  effects  (zero  rpm) .  These  pinned-pinned  modes  are 
used  in  later  calculations  to  form  the  coupled  modes. 


1  Numbers  in  brackets  refer  to  references  at  the  end  of  this  report. 
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Figure  4-3  Flow  Chart  Depicting  8660  Calculation  Procedure 
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Step  3  -  The  coupled  rigid  body  modes  are  computed  without  gyroscopics 
effects  by  assuming  the  substructures  to  be  Infinitely  stiff  and  joined  to 
each  ocher  through  the  linear  stiffness  elements.  Nonlinearities  are 
ignored  at  this  point  in  the  procedure. 

Step  4  -  The  coupled  system  modes  are  computed  without  gyroscopics  effects. 
The  previously  computed  substructure  uncoupled  plnned-pinned  modes  and  the 
coupled  rigid  body  modes  ara  used  as  degrees  of  freedom  in  this  calculation. 
The  resulting  coupled  system  modes  become  the  basis  for  all  subsequent  cal¬ 
culations. 

Step  5  -  Any  of  a  large  number  of  calculations  are  possible.  Among  those 
are  che  ones  associated  with  including  the  effects  of  squeeze  film  dampers. 
Specifically  they  are: 

e  Frequency  domain  response 

e  Time  domain  response 

•  Log  decrement  (damped  natural  frequency) 

The  damper  effects  are  included  in  the  nonlinear  response  calculations  as 
modal  forces.  They  are  included  in  the  damped  natural  frequency  calculation 
as  complex  stiffness  and  viscous  damping  coefficients.  This  treatment  is 
discussed  further  in  Section  V. 

This  discussion  was  presented  to  provide  a  general  description  of  the  rotor 
dynamic  analysis  computer  program  used  at  DDA.  The  models  are  developed 
using  a  lumped  mass,  beam  element  formulation.  Solutions  are  performed 
using  a  modal  synthesis  approach  and  nonlinear  effects,  such  as  supplied  by 
squeeze  film  dampers,  are  included  in  the  response  calculations  as  modal 
forces. 


3.  Damper  Component  Analysis  Tool 

The  damper  component  analysis  tool  used  in  forming  the  integrated  rotor/ 
damper  analysis  program  was  developed  by  MTI  and  was  delivered  to  DDA  for 
coupling  with  the  rotordynamics  code.  A  brief  description  of  this  tool  is 


presented  to:  1)  indicate  the  type  of  data  which  must  be  assimilated  in 
the  rotor  dynamics  program  and  2)  establish  some  of  the  equations  and 
terminology  used  in  the  integration  of  the  programs. 


General  Description 

The  squeeze  film  damper  analysis  is  centered  on  a  calculation  of  the  response 
of  the  film  to  an  imposed  velocity  at  a  particular  location.  A  finite 
difference  solution  of  the  Reynolds  equation  is  used  which  provides  algebraic 
difference  equations  linking  pressure  values  at  each  point  of  a  rectangular 
mesh  on  the  developed  damper  surface.  The  analysis  is  sufficiently  genera¬ 
lized  to  facilitate  both  circular  and  noncircular  orbit  calculations. 

The  squeeze  film  fluid  pressure  distribution,  in  its  most  general  form,  is 
dependent  on  the  bearing  journal  velocity  and  acceleration  components,  the 
squeeze  film  thickness,  cavitation  and  supply  pressure,  the  physical 
properties  of  the  film  fluid,  and  the  bearing  geometry.  The  Reynolds 
equation  for  the  squeeze  film  pressure  distribution  is  written  in  terms  of 
pressure  coefficients  of  velocity  and  acceleration.  This  provides  a  form 
of  the  component  pressures  which  are  nondimensionalized  and  tabulated  as  a 
function  of  journal  eccentricity.  Subsequent  damping  and  inertia  coefficient 
calculations,  which  are  functions  of  velocity  and  acceleration,  make  use  of 
the  tabulated  pressure  coefficients,  thereby  avoiding  the  need  to  re-solve 
the  Reynolds  equation. 

A  simplification  to  circular  orbit,  is  used  when  applicable,  to  further 
reduce  the  number  of  calculations  required  during  a  rotor  response  analysis. 
By  utilizing  the  relationship  between  velocity,  acceleration,  and  journal 
eccentricity,  a  table  of  damping  and  Inertia  coefficients  can  be  generated 
as  a  function  of  journal  eccentricity. 

The  comprehensive  squeeze  film  damper  analysis  accounts  for  the  following 
effects: 

•  Fluid  film  cavitation  is  treated  using  the  *  or  half-film 
approach  where  any  computed  subambient  pressures  are  set 
to  ambient,  thus  defining  the  cavitation  boundary. 


•  Fluid  Inertia  effects  are  accounted  for  by  adjusting  the 
Reynolds  solution.  These  effects  result  in  a  linear 
correction  to  the  viscous  damping  ti rm. 

e  Fluid  supply  conditions  are  determined  by  designating 
pressure  to  be  equal  to  inlet  pressure  over  the  region 
of  the  supply  hole  or  groove. 

•  Flow  is  determined  by  integrating  around  the  boundary  of 
interest.  The  local  flow  is  calculated,  using  the  com¬ 
puted  pressure  gradient,  viscosity,  and  film  thickness. 


•  The  damper  length  is  accounted  for. 


e  As  an  implicit  part  of  the  Reynolds  equation  solution, 
the  following  are  accounted  for: 

.  Finite  orbit  size  or  eccentricity 

.  Damper  clearance  and  diameter 

.  Fluid  viscous  effects. 


The  squeeze  film  damper  analysis  uses  a  fixed  Cartesian  x-y  coordinate  system 

with  origin  at  the  bearing  center,  as  shown  in  Figure  4-4.  In  this  system 

the  Journal  center  has  the  coordinates  x  and  y,  the  velocities  dx/dt  and 

2  2  2  2 

dy/dt,  and  the  accelerations  d  x/dt  and  d  y/dt  .  Points  in  the  lubricant 
film  are  located  by  the  axial  coordinate  z,  measured  from  the  centerplane,  and 
by  the  circumferential  angle  9,  measured  from  the  negative  x-axls. 


Figure  4-4  Squeeze  film  damper  model 


TE-80-123 


The  squeeze  film  damper  outer  cylinder  supports  a  non-rotating  bushing  on 
an  oil  film.  The  bearing  dimensions  and  lubricant  properties  are: 

R  ■  1/2D  ■  squeeze  film  inner  radius, (m*) 

L  ■  land  length, (m) 

C  ■  radial  film  clearance  when  centered, (m) 

2 

y  ■  lubricant  viscosity, (N. s/m  ) 

3 

p  »  lubricant  density , (kg/m  ) 
x,  y  ■  bushing  center  coordinates (m) 


0  *  circumferential  coordinate  angle,  measured  from 
negative  x-axis  (radians) 

z  »  axial  coordinate, (m) 

h  *  C  +  x  coscjt  +  y  sinuiT  *  local  film  thickness ,  (m) 

T  ■  time, (sec) 

3 

q^  ■  flow  from  k'th  feeding  hole,(m  /s) 
f  «  reference**  frequency , (Hz) 

2  «  2  »  f  *  reference**  angular  frequency, (rad /sec) 

2 

p  -  film  pressure, (N/n  ) 

Rotor  Response  Calculations 

Within  the  squeeze  film  analysis,  four  damping  coefficients  (B  ,  B  ,  B  , 

xx  xy  yx 

and  B  )  and  four  inertia  coefficients  (C  ,  C  ,  C  ,  and  C  )  are  deter- 
yy  xx  xy  yx  yy 

mined,  at  constant  supply  pressure,  for  a  given  damper  motion.  The 
instantaneous  forces  exerted  by  the  damper  during  orbit  are 


F  -  -B 

X  -B 

Y  -c 

•  • 

X 

-c 

•• 

Y 

X  XX 

xy 

XX 

xy 

e 

e 

M 

»# 

(4-1) 

F  -  -B 

X  -B 

Y  -C 

X 

-C 

Y 

y  yx 

yy 

yx 

yy 

*  The  units  shown  are  metric.'  Equivalent  English  units  also  apply. 

**  A  reference  frequency  is  used  eo  reduce  intermediate  pressure  tables  to 
non-dimensional  form. 
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Linear  Analysis  with  Circular  Orbit 
In  the  case  of  linear  circular  orbit 


B  -  B  ;  B 
yy  xx  yx 


-B  ;  C  -  C  ;  C 
xy  yy  xx  xy 


(4-2) 


These  damping  inertia  coefficients  are  used  directly  in  the  DDA  linearized 
rotor  dynamics  analysis  in  which  the  journal  orbit  is  treated  as  circular 
and  at  fixed  radius. 


Nonlinear  Analysis  with  Circular  Orbit 

The  circular  orbit  coefficient  tables  are  used  in  DDA  rotor  dynamics  analyses 
to  carry  out  nonlinear  steady-state  rotor  response  calculations  in  the  fre¬ 
quency  domain.  For  these  calculations  the  journal  orbit  radius  is  allowed 
to  vary  with  response  displacement,  thereby  introducing  variable  damping 
and  inertia  coefficients.  Thus,  the  instantaneous  forces,  as  defined  by 
equation  (4-1),  also  vary  with  response  displacement. 

Within  the  DDA  rotor  dynamics  analysis,  nonlinear  circular  orbit  calculations, 
using  the  squeeze  film  damper  analysis,  are  carried  out  iteratively.  An 
initial  response  value  of  radial  deflection  is  computed.  Using  this  deflec¬ 
tion,  the  damping  and  inertia  coefficients  are  determined.  These  coeffi¬ 
cients  are  used  to  form  the  real  space  forces  on  the  rotor  as  defined  by 
equation  (4-1) .  A  Newton-Raphson  technique  is  used  to  solve  for  the  next 
iteration  until  the  solution  converges.  This  is  fully  described  in  the 
Rotor/Damper  interface  section.  If  the  calculated  deflection  exceeds  the 
damper  clearance  the  unbalance  force  Is  reduced  and  then  increased  incre¬ 
mentally  with  convergence  at  each  step.  A  solution  is  attained  when  the 
unbalance  reaches  its  original  value. 


Nonlinear  Analysis  with  Noncircular  Orbit 

The  DDA  rotor  dynamics  analysis  uses  the  table  of  Reynolds  pressure  coeffi¬ 
cients  of  velocity  and  acceleration  when  the  nature  of  the  problem  demands 
that  noncircular  orbits  be  considered.  In  cases  such  as  the  time  history  of 
nonsynchronous  response  to  unbalance,  the  journal  orbit  is  noncircular  and 
the  journal  center  velocity  and  acceleration  components  are  not  directly  a 


86 


function  of  radial  displacement.  The  damping  and  inertia  coefficients  are 
calculated  based  on  pressure  coefficients  interpolated  using  the  Reynolds 
pressure  coefficient  table.  Thus,  the  forces  defined  in  equation  (4-1)  are 
calculated  for  each  rotor  state  ($,  fc,  Y ,  X,  Y)  in  a  given  time  step. 

Within  each  time  step  these  forces  are  equilibrated,  and  the  orbital 
trajectory  is  determined  using  a  predictor-corrector  technique. 


SECTION  V 


ROTOR/DAMPER  INTERFACE 

Squeeze  film  dampers  interface  with  the  rotor  dynamic, analysis  program  in 
the  following  calculations: 

e  Linear 

-  Undamped  critical  speeds  -  (Define  need  for  squeeze  film 

damping) 

-  Log  decrement  -  (Define  squeeze  film  damping  required  for 

maximum  modal  damping) 

-  Frequency  domain  response  -  (Define  required  squeeze  film 

damper  motion) 

•  Nonlinear 

-  Frequency  domain  response  -  (Define  overall  unbalance 

response) 

-  Time  domain  response  -  (Define  suddenly  applied  load 

effects  and  2  plane  non-circular  steady  state 
response) 

The  linear  calculations  are  performed  using  circular  orbit,  linearized 
coefficients.  The  non-linear  calculations  are  performed  using  either 
circular  orbit  forces,  computed  from  the  coefficients,  or  non-circular 
orbit  forces,  computed  from  the  velocity  and  acceleration  pressure  distri¬ 
bution  tables.  Each  of  these  areas  of  calculation  is  discussed. 

1.  Linear  Calculations 

Squeeze  film  dampers  are  inherently  non-linear  devices  as  can  be  seen  in  the 
damper  force  plot  of  Figure  5-1.  This  characteristic  was  computed  for 
circular  whirl  of  a  4.5-inch  diameter  bearing  with  a  length-to-diameter  (L/D) 
ratio  of  0.11  and  a  radial  clearance  C  of  0.005  inches.  The  tangential 
(damping)  force  is  relatively  linear  for  an  eccentricity  ratio  (e/C), 
between  0.0  and  0.6.  For  eccentricity  ratios  larger  than  0.6  the  damping 
increases  rapidly.  The  radial  (stiffness)  force  component  is  substantially 
lower  than  the  tangential  force  through  an  eccentricity  ratio  of  about  0.5. 
However,  it  is  linear  only  through  about  0.4.  The  radial  force  also  has 
a  stiffening  characteristic  for  higher  eccentricity  ratios.  Although  the 
damper  is  clearly  non-linear,  certain  information  related  to  critical  speeds, 
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The  damper  component  computer  program,  SQZDMP,  described  in  earlier  sections 

of  this  report,  can  produce  stiffness  and  damping  Information  in  several 

forms.  Circular  orbit  coefficients,  computed  as  a  ratio  of  tangential  and 

radial  force  to  tangential  velocity  for  purely  circular  whirl  are  presented. 

These  are  the  direct  damping  (B  )  and  cross-coupling  damping  (B  )  coeffl- 

xx  xy 

dents.  They  are  used  in  the  following  equation  for  the  tangential  and 
radial  forces. 


where 


F  -  tangential  force  (damping) 

-  radial  force  (stiffness) 
e  -  radial  eccentricity  (deflection) 
i 'i  -  whirl  velocity 


If  an  eccentricity,  e  and  whirl  speed,  are  assumed,  then  for  circular  whirl 

(5-2) 


F  «  B  e  't>  -  (B  e)  '4 

t  XX  XX 


F  *  B  e  2  *  (B_e)  <6 
r  xy  xy 


or 


-  (1  +  ig)K 


K  -  B  55 
xy 

a  m  B  /B 
*  XX  xy 


where 

Kp  -  dynamic  stiffness  coefficient  (lb/in.) 
K  -  spring  stiffness  coefficient  (lb/in.) 
g  -  structural  damping  coefficient 


(5-3) 


This  allows  the  treatment  of  a  damper  as  a  dynamic  stiffness  element  with 
linear  radial  stiffness  and  structural  damping  for  an  assumed  circular 
orbit  eccentricity  and  speed.  This  simplifies  the  numercial  calculations 
and  is  only  valid  for  linear  steady  state  response  calculations. 


Another  approach  treats  the  damper  as  both  a  stiffness  and  a  viscous  damping 
element.  Again,  assuming  a  circular  orbit  eccentricity  of  e  and  a  whirl 
speed  of  *4. 
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or 


K  -  B 
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where 

K  -  spring  stiffness  coefficient  (lb/in.) 

C  -  viscous  damping  coefficient  (lb  sec/ln.) 


(5-5) 


When  the  whirling  motion  can  no  longer  be  considered  as  a  centered,  circular 
orbit,  the  problem  is  more  complicated.  However,  there  still  is  a  need  to 
evaluate  the  damper  on  a  basis  much  simpler  than  becoming  Involved  in  a  time 
history  analysis.  Therefore,  provisions  have  been  made  in  the  damper 
component  analysis  to  compute  effective  stiffness  and  damping  values  for 
non-circular  orbits.  These  have  been  computed  from  a  hysteresis  loop  rela¬ 
tionship  between  force  and  displacement  in  orthogonal  directions  as  indicated 
in  Figure  5-2.  The  effective  linear  stiffness  in  the  X  direction,  K^,  is 
the  slope  of  the  line  joining  points  of  maximum  and  minimum  displacement. 

The  effective  linear  damping  coefficient  is  determined  from  the  area 
enclosed  by  the  hysteresis  loop  and  the  relationship  indicated.  This  will 
produce  a  viscous  damping  coefficient,  B^.  Similar  effective  linear 
coefficients  are  computed  for  the  T  direction,  and  B^.  These  may  be 
used  in  an  asymmetric  linear  calculation  or  averaged  for  an  assumed  circular 
analysis.  Structural  damping  coefficients  may  be  obtained  from  the  viscous 
coefficients  by  the  relationship. 

*x  -  Bx 

(5-6) 

*y  -  Bx  */Ky 

where  if  is  the  whirl  speed  in  rad/sec. 

Each  of  these  approaches  to  determining  damper  stiffness  and  damping  coeffi¬ 
cients  is  assimilated  by  computer  program  BB60  in  open  loop  fashion  (l.e. 
there  is  no  computerized  Interface  with  these  calculations).  However,  the 
ability  to  evaluate  damping  elements  with  these  characteristics  exists  and 
are  briefly  described. 
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Undamped  Critical  Speeds 

As  previously  stated  in  Section  4.0,  all  rotordynamics  calculations  are 
performed  using  a  set  of  coupled  system  modes  at  zero  speed  as  the  basic 
degrees  of  freedom.  These  modes  were  generated  from  the  modal  synthesis 
of  substructure  pinned-pinned  modes  and  coupled  rigid  body  modes.  The 
rigid  body  modes  consider  the  systems  to  be  rigid  and  coupled  by  inter¬ 
connecting  stiffnesses.  A  dynamic  matrix  is  formed  with  gyroscopic  effects 
included  in  the  mass  matrix.  This  is  done  by  setting  a  speed  ratio.  This 
is  the  ratio  of  the  undamped  mode  frequency  to  the  system  rotor  speed.  The 
degrees  of  freedom  are  the  linear  complete  engine  modes.  The  critical 
speeds  and  mode  shapes  are  calculated  for  the  rotor  operating  at  specified 
ratios  to  the  critical  speeds.  For  a  single-rotor  engine,  the  speed  ratio 
would  be  1.0.  A  two-rotor  system  would  have  two  sets  of  critical  speeds. 

For  example,  an  engine  with  a  counterrotating  rotor  with  rated  rotor  speeds 
of  10,000  and  7,500  rpm  would  have  speed  ratios  sets  of  (1.0,  -0.75)  and 
(-1.33,  1.0).  The  first  set  of  modes  are  excited  by  the  first  rotor,  and 
the  second  set  of  modes  are  excited  by  the  second  rotor.  The  analysis  Is 
performed  in  modal  space,  and  the  results  are  converted  to  real  space  before 
output.  Frequency-speed  maps  can  be  generated  with  this  calculation  pro¬ 
cedure  by  performing  the  calculation  at  various  speed  ratios.  The  calcula¬ 
tion  will  allow  the  estimate  of  damper  effects  on  critical  speeds  and  mode 
shapes  as  a  first  order  evaluation. 

Logarithmic  Decrement 

The  logarithmic  decrement  (log  dec)  calculation  is  used  to  determine  the 
mount  damping  required  for  maximum  modal  damping.  The  damped  modal  natural 
frequencies  and  associated  log  dec  may  be  plotted  as  a  root  locus  with  mount 
damping  varied.  This  information  is  used  to  select  the  amount  of  desirable 
squeeze  film  damping. 

This  calculation  is  specifically  used  to  determine  the  effect  of  direct  and 
indirect  damping  on  mode  damping  and  frequencies.  Two  options  are  available 
for  log  dec  calculations.  The  first  accepts  a  direct  viscous  damping  input  and 
calculates  log  dec.  The  second  accepts  a  damper  eccentricity  and  interpolates 


93 


for  both  direct  and  indirect  damping.  Input  of  damper  eccentricity  is 
included  on  the  squeeze  film  damper  input  card  to  BB60.  The  method  for 
calculation  is  the  same  in  both  cases  and  is  given  below. 

The  log  dec  calculation  involves  the  solution  of  the  following  modal 
differential  equations. 

[m]  <h)  +  [c]  <h)  +  0]  (6)  +  [lG]  (h)  -  (o)  (5-7) 

where 

Qf]  “  modal  mass  matrix 

(Y)  “  complex  stiffness  matrix  substructure  with  mount 
structural  damping 

[c]  ■  viscous  damping  matrix  in  modal  degrees  of  freedom 
flGJ  ■  gyroscopic  inertia  matrix 
Q  ■  rotor  speed 

(h)  •  vector  of  modal  participations 

Each  modal  degree  of  freedom  has  associated  with  it  a  normalized  deflected 
mode  shape.  The  coupled  modes  will  Include  participation  of  all  uncoupled 
modes.  The  complex  vector  (h)  will  be  the  participation  vector  of  the 
uncoupled  modes. 

The  solution  of  the  coupled  modal  differential  equation  is  obtained  by 
letting  h  ■  hoe*C. 

where 

X  -  complex  eigenvalues  ■  r  +  iw 
r  -  modal  damping 
w  -  damped  natural  frequency 

The  characteristic  equation  becomes: 

+  »H<0 +  0c3 +  S3|- »  '5-<» 


1 


3 
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The  energy  dissipation  for  the  mounts  will  be  transformed  from  local  to 
modal  degrees  of  freedom  to  be  Included  In  the  modal  damping  matrix 


GO  •  BtCcsfD  H] 

where  fC---]  -  squeeze  film  damping  matrix  in  local  degrees 
of  freedom  (X) 

X  *  (jfJCh) ,  QTj  -  is  the  transformation  matrix 


(5-9) 


The  log  dec  will  be  calculated  for  each  mode  as:  6  ■  -  —  where  w  ^  0 

V 

is  in  RAD/ SEC. 


Steady  State  Response 

The  results  of  this  analysis  are  the  forced  deflections  at  each  of  the  mass 
station  due  to  an  imposed  unbalance  distribution.  With  the  optimum  mount 
damping  input,  the  resulting  response  at  the  mount  is  that  expected  for  the 
damper  or  dampers.  To  keep  the  damping  in  the  linear  range  the  target 
squeeze  film  radial  clearance  is  set  to  1.7  times  the  1/2  peak-to-peak 
response . 

The  input  damping  is  included  as  part  of  the  viscous  damping  matrix  or  as  an 
imaginary  stiffness  component.  For  steady  state  response  the  damping  may  be 
included  using  either  approach.  In  BB60  it  is  included  as  an  imaginary 
stiffness  ratio  (g)  or  as  a  viscous  damper  (c) . 

The  equation  solved  is  written  as  follows: 

[M]  (q)  +  [C]  (q)  +  [X]  (q)  -  (f)  (5-10) 

ft 

where  M  is  the  mass  matrix  with  gyroscopic  terms  included  as  previously 
discussed. 

Q]J  is  the  damping  matrix 

QQ  is  the  complex  stiffness  matrix 

( q is  the  displacement  degrees  of  freedom 

{ is  the  forcing  function  due  to  rotating  unbalance 


The  zero  frequency  undamped  mode  shapes  are  calculated  using  the  real  part 
of  the  stiffness  matrix  (k^J  and  the  mass  matrix  without  gyroscopic  terms. 
The  transformation  from  displacement  to  undamped  zero  frequency  modal  degrees 
of  freedom  is: 

(?)  -  CO  0}  <5-n> 

{h^  -  will  be  referred  to  as  the  modal  displacement  participation  fac- 
tors  in  the  remaining  discussion. 


The  dynamic  matrix  may  be  written  as  follows: 

(?D  5  H  +  n0  1  "  n2  0]  (5-12) 

where  ft  is  the  rotor  speed  at  which  the  unbalance  occurs. 


CO  (?)  -  0)  can  be  written  for  modal  degrees  of  freedom  as: 

ETGD&3 00 •  07  (») 


(5-13) 


This  reduces  the  computational  complexity. 


2.  Nonlinear  Calculations 

The  linear  calculations  just  described  are  useful  in  screening  damper  designs 
and  allow  the  elimination  of  damper  variables  and  configurations  which  are 
unsatisfactory.  It  provides  an  analysis  which  can  be  used  to  "get  in  the 
neighborhood"  of  a  solution.  Nonlinear  analyses  must  then  be  conducted  to 
further  quantify  the  effects  of  squeeze  film  dampers  on  the  rotor  dynamic 
behavior  of  a  rotor  system.  The  damper  component  software,  described  in 
earlier  sections,  has  been  combined  with  the  rotor  dynamic  software, 
described  in  Section  IV,  to  provide  coupled  evaluations  of  dampers  in  gas 
turbine  engines.  The  resulting  coupled  response  analyses  In  the  frequency 
and  time  domain  are  discussed. 

Steady  State  Response 

This  routine  calculates  steady  state,  circular  orbit  response  for  specific 
squeeze  film  designs  and  is  used  to  make  damper  refinements. 
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The  coupled  rotor /damper  analysis  for  frequency  domain  response  is  per¬ 
formed  with  the  squeeze  film  damper  forces  Included  as  external  forcing 
functions  which  are  dependent  on  mount  displacement.  Therefore,  the 
solution  is  iterative.  The  squeeze  film  damper  forces  are  a  function  of 
absolute  displacement  and  are  resolved  into  radial  and  tangential  components 
for  use  in  BB60.  The  complex  dynamic  matrix  that  was  calculated  for  the 
linear  response  problem  is  now  expanded  for  real  and  Imaginary  response 
components.  The  reason  for  this  will  be  shown  later. 

[■»]  0)  ■  DO*  DO  CO  0}  ■  Kf  (0  -  (?) 

[»■>]  •  [“0  +  Q»  J  1 

*  (5-14 

{"}  ■  (V  +  (hi) 1 

(?)  •  (Or)  +  (*i)  1 


IDDr]  \)  -  CDDl]  (hi)  ■  («r) 

C“l  J  (\)  +  ?»r]  (hl)  -  («l) 


(5-15) 


The  generalized  force  vectors  for  the  nonlinear  response  calculation  Includes 
the  squeeze  film  reaction  forces  f  ( jVj)  •  These  forces  are  stored  in  tables 
of  the  radial  and  tangential  forces  versus  displacement  eccentricity  (£qj / C) , 
where  C  is  the  squeeze  film  clearance.  The  radial  f  ( [q  )  and  the  tangential 

K  w  ~ 

^T  ^  "Vj)  ^orces  are  rotated  into  BB60's  real  displacement  coordinates  by  the 
following  equations:  F({q])  -  fR(QT]>  +  ifT((jfj> 

Ke(B([q]))  -  fR([q])  ^  -  VK)>  (5'16) 

^  W  J 

and 

I»(B(0|]))  -  fR([qj  ^ 

«)■  (9r)- {«-<&]))  (5-17) 

(qi)‘(Qi)-  CpID(Im»<C’])) 
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(qj)  and  (qJ)  replace  (Q,)  «ad  (Qj) 
on  right  Hand  aid*  of  equations. 

Sine*  the  squeeze  film  force*  *r*  iteratively  determined  from  the  relative 
displacement  of  the  journal,  the  equation  will  not  exactly  balance  and  a 
residual  force  ( 00)}  may  be  calculated 

(fs  'S’)*-  Cdd*3  i\)  -  M  (hi) +  &TD  {«*»'&>;  -  (V 


(5-18) 


(fA  ([q])}  x  -  Q>Dj]  (hR)  +  [ddJ  (hT)  +  [pT]  {lmB([q])}  -  (Qj. 

To  converge  the  residual  forces  to  zero  a  Newton  Raphson  iteration  is  used 

- ~  llyk.ii;  •  (W -  ("J + 


ft  ( nr ))  •  (~dd 

(.VLilV  ^fA((q|)I)  L- 


(Kit)} 

a(q) 


where 

[j<B<|ql)>] 


} 


[dd]t  +  [pj  jj(B(|q{))]  [Pj 


>r] 

—  i 

M  i 

1  1 
*-» 

ii 

1 

(5-19) 


Sum  over  all  non¬ 
linear  mounts 


3ReB(q)  „  _ 1 

^  kl 


~3ReB 

3ReB 

3  Req 

3  Imq 

3Imb  ' 

3lfflB 

_3Req 

3  Imq  - 

(Rea)^ 

2l3|q| 

(5-20) 


Req  -  Real  part  of  q 
Imq  -  Imaginary  part  of  q 


.[hi  - 

c*00*)  [ftl  -  Sf]} 


(5-21) 
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(5-22) 


3ReB(q) 

3Imq 


3ImB(q) 

3Req 


aimB(q)  _ 
3Imq 


J 


(5-23) 


(5-24) 


The  terms  do  not  satisfy  the  Cauchy-Riemann 

functions.  This  means  that  3ReB(  [gQ  i 

3Req 


conditions  for  analytic  complex 

3Im(B(  |q|)) 

3Imq 


and 


(5-25) 


aReBliSU  a 

3Imq 


-3ImB( |g I ) 
3Req 


Consequently,  the  matrix 


was  expanded  into  real  and  imaginary  parts. 


The  iteration  is  begun  by  assuming  the  squeeze  film  forces  to  be  zero  and 

calculating  the  displacement  of  the  mount.  If  the  displacement  is  less  than 

dB  dB 

the  damper  clearance  the  derivatives  xx  and  xx  are  calculated  from  a 

de  de 

spline  fit  of  the  tabular  data  in  the  vicinity  of  the  calculated  displace¬ 
ment  'q"'.  The  eccentricity  e  is  equal  to  Tq^./C. 

<w  U  * 

If  the  damper  clearance  c  is  exceeded,  then  the  unbalance  forces  are 
reduced  and  then  incremented  upward  to  reach  the  original  unbalance.  The 
final  values  of  the  state  variables  are  used  as  starting  values  for  the  next 
iteration. 
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The  radial  and  tangential  forces  and  used  in  the  iteration  are  cal¬ 
culated  as  follows: 

2 

f _  »  -  Bxxft  q  +  Cxyfl  q 

(5-26) 

2 

fT  -  -  Bxyfl  q  +  Cyyfl  q 

This  analysis  assumed  circular  orbits  so  that  f  and  f  are  steady  rotating 

K  X 

forces. 


The  radial  and  tangential  force  derivatives  are  given  below: 


dB 

_ 

diet 


They  are  written  with  the  clearance  (c)  factored  out. 


(5-27) 


The  coordinate  system  used  for  the  rotor  dynamic  analysis  is  left-handed  and 
is  right-handed  for  the  damper  component  analysis.  Figure  5-3  compares  the 
two  coordinate  systems.  Figure  5-4  gives  an  overview  of  the  non-linear 
steady  state  calculation. 


Time  Domain  Response 

So  far  the  discussion  has  covered  the  case  of  coupled  rotor/damper  analyses 
having  orbits  with  circular  symmetry.  The  general  case  of  damper  motion  is 
noncircular  and  non-centered  in  the  damper  clearance  space.  The  proper 
approach  to  evaluating  this  condition  is  to  use  a  time  domain  solution  so 
that  the  non-circular  orbit  can  be  tracked  as  a  function  of  time  and 
noncircular  state  variables  and  damper  forces  can  be  computed.  DDA  has 
modified  it's  existing  rotor  dynamic  transient  analysis  to  incorporate  the 
effects  of  a  squeeze  damper  for  both  circular  and  non-circular  whirl.  The 
following  is  a  discussion  of  the  DDA  time  domain  response  calculation 
followed  by  discussions  of  the  circular  orbit  and  non-circular  orbit  inter¬ 
faces. 
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c)  Coordinate  Transformations 


Figure  5-3.  Comparison  of  Component  and  System  Coordinate  Systems. 
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ENTER 

PROGRAM 


BB60  INPUT  DATA 


Figure  5-4  Overview  of  non-linear  steady  state  response  calculation 
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Response  Calculation 

This  analysis  calculates  the  response  of  the  engine  to  various  shock 
excitation  such  as 

•  Initial  velocity  (a  drop) 

•  Half  sine  displacement 

•  Rectangular  displacement 

•  Series  displacement 

•  Verslne  displacement 

•  Half  sine  acceleration 

•  Rectangular  acceleration 

•  Series  acceleration 

•  Versine  acceleration 

•  Suddenly  applied  unbalance 

The  displacement  shocks  are  displacements  of  the  ground  while  the  accelera¬ 
tion  shocks  are  applied  to  the  engine.  The  series  shock  input  is  a  fourth 
order  polynomlnal.  The  option  of  most  initial  interest  for  rotor  dynamics 
is  the  time  response  to  suddenly  applied  unbalances.  It  is  with  this  option 
that  evaluations  will  be  made  as  to  the  coupled  effect  of  squeeze  film 
dampers  on  rotor  dynamic  behavior  during  conditions  of  blade  loss.  If  the 
transient  solution  is  allowed  to  decay,  the  steady  state  solution  results. 
This  method,  which  considers  the  effects  of  gravity,  is  a  more  accurate 
model  of  rotor/damper  systems  which  have  offset  noncircular  orbits.  An 
additional  feature  also  allows  the  evaluation  of  the  rotor /damper  system 
using  circular  orbit  coefficients.  The  solution  is  capable  of  analyzing 
either  linear  or  nonlinear  systems. 

The  calculation  is  performed  in  modal  space  with  a  time  marching  integration. 
A  predictor-corrector  method  with  a  variable  time  step  is  used.  The  basic 
equations  are: 
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J 


] 

> 


'4 


(*)  •  (X)l-1  +  «  C1'5®!-!  - 

(x)  ■  (x)i-i +  -5  “  C(x)i-i +  (x)i] 


Predictor  ■  [V}  ^  (f^ 


(f)  -  (fl)  cos  jSj  +  (Fv)  sin  ^  -  (x)(x)  +  (ij  -  (3(x}  ♦  (l*) 

(*}i  ■  (X)i-1  +  4t  [?  (*}l-l  '  if  (X)i-2  +  if  CS)i]  "I 


(5-28) 


(X)i  “  (X)i-1  +  At  [f  (X)i-1  *  I?  (X)i-2  +  16  COil  J 

*—  -J 


Corrector 


where 


(F) 

(fl) 

(Fv) 

00 

B 

(F.) 

CO 

(L«) 


00 


is  the  modal  force  vector 

is  the  modal  lateral  force  due  to  unbalance 

is  the  modal  vertical  force  due  to  unbalance 

is  the  modal  displacement 

is  the  diagonal  modal  stiffness  matrix 

is  the  modal  force  due  to  acceleration 

is  the  damping  gyroscopic  matrix 

is  the  modal  forces  due  to  nonlinearities 

This  is  calculated  at  each  step  from  the  displacements  and 
nonlinear  mounts. 

is  the  diagonal  modal  mass  matrix 


The  Interface  between  this  calculation  in  the  rotor  dynamic  program  and  the 
damper  component  program  is  made  through  the  modal  force  vector  L^.  The 
damper  reactions  at  each  increment  of  time  are  treated  as  modal  forces.  The 
modal  forces  need  to  be  re-evaluated  at  each  of  these  increments  in  terms  of 
the  real  space  state  variables  (i.e.  vertical  and  lateral  displacement, 
velocity,  and  acceleration).  This  is  described  in  the  section  called 
Non-Linear  Steady  State  Response  Calculation  in  some  detail  along  with  the 
dynamic  matrix  equation.  The  state  variables  are  passed  to  interpolation 
routines,  one  for  circular  orbit  (1  plane)  and  another  for  noncircular  orbit 
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(2  plane) .  Again  the  calling  state  variables  are  transformed  to  account 
for  differences  in  the  coordinate  system  as  seen  in  Figure  5-3.  The  state 
variables  are  defined  as: 


Damper 

Variables 

*D 

*D 

*D 


Rotor 

Variables 


With  this  transformation  the  forces  returned  from  the  interpolation  routines 
for  rotor  program  use  are  -F^  an<*  As  in  tlle  circular  orbit,  nonlinear, 

steady  state  calculation,  the  real  space  forces  are  transformed  to  modal 
forces  1^  by  the  equation 

{*»>•  COfr.)  (5-29> 

where,  as  before 

matrix  of  mode  shapes 
real  space  complex  force  vector 

This  process  is  repeated  at  each  time  step  in  the  simulation. 


Circular  Orbit  Interface 

The  Interface  with  the  circular  orbit  option  of  the  time  domain  response 
calculation  is  made  through  a  call  to  the  interpolation  subroutine  FIND  in 
the  same  manner  described  for  the  circular  orbit,  steady  state  analysis. 
rTND  Interpolates  the  circular  orbit,  normalized  coefficient  tables.  The 
coefficients  are  converted  to  real  dimensions  and  the  F^  and  F^  forces 
are  computed  using  equation  (5-29).  As  before,  the  forces  which  are  returned 

and  Yr 

orthogonal  axes  and  are  respectively,  -F,_  and  F  _.  The  only  difference 


to  the  calling  routine  in  BB60  are  opposing  motion  along  the 


between  this  interface  and  the  Interface  with  the  steady  state,  circular 
orbit  calculation  is  that  nonzero  values  of  the  rotor  velocity  and 
acceleration  are  used.  The  important  aspect  of  this  interface  relates  to 
the  difference  in  the  rotor  dynamic  and  damper  component  coordinate  systems. 
As  a  result,  two  coordinate  transformations  are  required.  The  rotor  dynamic 
state  variables  must  be  transformed  to  the  damper  component  coordinate 
system  prior  to  calling  subroutine  FIND.  This  permits  the  calculations  to 
occur  in  the  damper  analysis  in  its  programmed  coordinates.  However,  upon 
the  return  of  forces  to  the  rotor  dynamic  analysis,  the  forces  must  be 
transformed  from  the  damper  coordinate  system  to  that  of  the  rotor.  This 
procedure  is  indicated  in  Figure  5-4.  Figure  5-5  gives  an  overview  of  the 
non-linear  time  transient  analysis. 

Noncircular  Orbit  Interface 

The  interface  with  the  noncircular  orbit  option  of  the  time  domain  response 
calculation  is  made  through  a  call  to  the  interpolation  subroutine  FORCE. 
Values  called  to  the  subroutine  from  the  damper  analysis  are  values  of  the 
state  variables  which  have  been  transformed  to  account  for  differences  in 
the  damper  and  rotor  coordinate  systems  (Figure  5-3).  Within  the  sub¬ 
routine  the  damper  forces  F^  and  F^  of  the  damper  component  analysis  are 
computed  from  stored  pressure  data.  These  data  assume  the  bearing  to  have 
axlsymmetrlc  geometry  so  that  the  only  parameter  is  eccentricity  e.  The 
stored  data  refer  to  a  rotating  coordinate  system  and  are  used  to  compute 
the  rotating  damper  forces.  These  are  then  transformed  to  the  fixed  X^— 
coordinate  system.  The  values  of  F^  and  F^  are  then  returned  to  the 
calling  routine  as  -F,^  and  F^R  respectively  for  BB60  use. 


EXIT  PROGRAM 


Figure  5-5  Overview  of  non-linear  time  transient  analysis 


SECTION  VI 


SAMPLE  ROTORDYNAMI CS  PROBLEM 

The  analytical  tools  that  have  been  described  are  demonstrated  here  for  a 
typical  jet  engine.  A  system  geometric  representation  Is  given  in  Figure 
6-1  which  shows  the  O.D.  and  1.0.  plotted  versus  engine  station  for  the 
rotor  and  case.  This  information  is  used  to  develop  the  engine  stiffness 
model.  Masses  and  inertias  are  lumped  at  locations  along  the  rotor  to 
represent  the  compressor  and  turbine  stages.  Masses  are  also  distributed 
along  the  case.  The  rotor  is  mounted  to  the  case  at  engine  stations  0.0  and 
65.70.  The  case  is  also  mounted  to  ground  at  stations  0.0  and  65.70.  The 
aft  rotor  mount  has  a  150,000  lb/in.  isolator  in  series  with  the  bearing. 

The  forward  mount  is  in  series  with  a  squeeze  film  damper.  A  100,000  lb /in. 
isolator  is  in  parallel  with  the  damper.  The  damper  film  is  10.34  inches 
in  diameter  and  its  width  is  2.59  inches.  The  oil  is  fed  through  unsealed 
ends. 

The  sample  analyses  shown  here  include  calculations  of: 
e  Critical  Speeds 
e  Log  Dec 

•  Response  to  Unbalance  (Steady  State) 

-  Normal 

-  Abusive 

•  Time  Transient  (Shock  Loading) 

-  Small  Unbalance  (Normal) 

-  Blade  Loss  (Abusive) 

The  critical  speeds  and  mode  shapes  for  this  sample  problem  are  shown  in 
Figure  6-2  for  a  front  rotor  to  case  stiffnesses  of  500,000  and  100,000  lb/in. 
These  plots  may  be  made  automatically  when  requested  in  the  rotor  dynamic 
input . 

Figure  6-3  shows  the  damped  eigenvalues  for  the  typical  jet  engine  with 
front  compressor  mount  viscous  damping  varied.  The  parallel  stiffness  is 
held  constant  at  100,000  lbs/ln.  This  shows  the  effect  of  direct  damping 
only. 
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ROTOR  (COMPRESSOR  AND  TURBINE) 

WEIGHT  -  810  LBS 
C.G.  -  37.4  IN 

PITCHING  INERTIA  ABOUT  C.G.  -  1172  LB-$EC2/IN 
MOUNTED  FRONT  -  100,000  LBS/IN  ISOLATOR  +  SQUEEZE  FILM  DAMPER 
BACK  -  150,000  LBS/IN  ISOLATOR 
1ST  PINNED  MODE  -  9,900  CPM 

CASE 

WEIGHT  -  1687  LBS 
C.G.  -  29.8  IN 

PITCHING  INERTIA  ABOUT  C.G.  -  3528  LB-$EC2/IN 
MOUNTS  FRONT  -  200,000  LBS/IN 
BACK  -  50,000  LBS/IN 
1ST  PINNED  MODE  «  7,370  CPM 


Figure  6-1  Sample  problem  input  description 


Figure  6-3 

Sample  Problem  for  Log  Dec  Calculation,  Direct  Damping  Varleo 
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Figure  6-4  shows  the  damped  eigenvalues  with  both  direct  and  indirect 
viscous  damping  varied.  The  variable  is  the  squeeze  film  damper  eccentricity. 

The  log  dec  calculation  just  demonstrated  would  be  used  to  help  size  and 
understand  the  operation  o£  the  squeeze  film  damper. 

The  steady  state  response  is  shown  in  Figures  6-5  and  6-6  for  a  normal 
unbalance  of  2  ounce  inches  and  for  a  blade  loss  with  19.5  in/lbs  unbalance 
2.28  inches  aft  of  the  damper.  Automatic  plotting  of  response  vs.  speed  is 
available  if  requested.  The  rotor  time  transient  response,  for  abusive 
unbalance,  at  the  squeeze  film  damper  is  shown  in  Figure  6-7  for  a  damper 
with  a  clearance  of  7  mils.  The  transient  analysis  uses  circular  orbit 
coefficients.  The  case  response  at  the  squeeze  film  is  also  given.  The 
case  is  responding  mostly  at  its  two  lowest  rigid  body  frequencies  of  20  Hz 
and  40  Hz.  Figure  6-8  shows  the  transient  response  for  a  normal  unbalance. 
These  analyses  just  demonstrated  would  be  typically  used  to  design  a  squeeze 
film  damper. 


Figure  6-5  Steady  State  Response  vs.  Station 
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Figure  6-7  Time  Transient  Response  to  Abusive  Unbalance 
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